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ABSTRACT 

We study the implications of Planck data for models of dark energy (DE) and modified gravity (MG), beyond the standard cosmological constant 
scenario. We start with cases where the DE only directly affects the background evolution, considering Taylor expansions of the equation of 
state w(a), as well as principal component analysis and parameterizations related to the potential of a minimally coupled DE scalar field. When 
estimating the density of DE at early times, we significantly improve present constraints and find that it has to be below « 2 % (at 95% confidence) 
of the critical density even when forced to play a role for z < 50 only. We then move to general parameterizations of the DE or MG perturbations that 
encompass both effective field theories and the phenomenology of gravitational potentials in MG models. Lastly, we test a range of specific models, 
such as k-essence, f(R) theories and coupled DE. In addition to the latest Planck data, for our main analyses we use background constraints from 
baryonic acoustic oscillations, type-la supemovae and local measurements of the Hubble constant. We further show the impact of measurements 
of the cosmological perturbations, such as redshift-space distortions and weak gravitational lensing. These additional probes are important tools 
for testing MG models and for breaking degeneracies that are still present in the combination of Planck and background data sets. 

All results that include only background parameterizations (expansion of the equation of state, early DE, general potentials in minimally-coupled 
scalar fields or principal component analysis) are in agreement with ACDM. When testing models that also change perturbations (even when 
the background is fixed to ACDM), some tensions appear in a few scenarios: the maximum one found is ~ 2 tr for Planck TT+lowP when 
parameterizing observables related to the gravitational potentials with a chosen time dependence; the tension increases to at most 3 cr when 
external data sets are included. It however disappears when including CMB lensing. 
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1. Introduction 

The cosmic microwave background (CMB) is a key probe of 
our cosmological model (Planck Collaboration XIII 2015), pro¬ 
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viding information on the primordial Universe and its physics, 
including inflationary models (Planck Collaboration XX 
2015) and constraints on primordial non-Gaussianities 
(Planck Collaboration XVII 2015). In this paper we use 


1 




Planck Collaboration: Planck 2015 results. XIV. Dark energy and modified gravity 


the 2015 data release from Planck 1 (Planck Collaboration I 
2015) to perform a systematic analysis of a large set of dark 
energy and modified gravity theories. 

Observations have long shown that only a small fraction of 
the total energy density in the Universe (around 5%) is in the 
form of baryonic matter, with the dark matter needed for struc¬ 
ture formation accounting for about another 26 %. In one sce¬ 
nario the dominant component, generic ally referred to as dark 
energy (hereafter DE), brings the total close to the critical den¬ 
sity and is responsible for the recent phase of accelerated ex¬ 
pansion. In another scenario the accelerated expansion arises, 
partly or fully, due to a modification of gravity on cosmologi¬ 
cal scales. Elucidating the nature of this DE and testing General 
Relativity (GR) on cosmological scales are major challenges 
for contemporary cosmology, both on the theoretical and ex¬ 
perimental sides (e.g., LSST Science Collaboration et al. 2009; 
Amendola et al. 2012a; Clifton et al. 2012; Joyce et al. 2014; 
Huterer et al. 2015). 

In preparation for future experimental investigations of DE 
and modified gravity (hereafter MG), it is important to determine 
what we already know about these models at different epochs 
in redshift and different length scales. CMB anisotropies fix the 
cosmology at early times, while additional cosmological data 
sets further constrain on how DE or MG evolve at lower red- 
shifts. The aim of this paper is to investigate models for dark 
energy and modified gravity using Planck data in combination 
with other data sets. 

The simplest model for DE is a cosmological constant. A, 
first introduced by Einstein (1917) in order to keep the Universe 
static, but soon dismissed when the Universe was found to be ex¬ 
panding (Lemaitre 1927; Hubble 1929). This constant has been 
reintroduced several times over the years in attempts to explain 
several astrophysical phenomena, including most recently the 
flat spatial geometry implied by the CMB and supernova obser¬ 
vations of a recent phase of accelerated expansion (Riess et al. 
1998; Perlmutter et al. 1999). A cosmological constant is de¬ 
scribed by a single parameter, the inclusion of which brings the 
model (ACDM) into excellent agreement with the data. ACDM 
still represents a good fit to a wide range of observations, more 
than 20 years after it was introduced. Nonetheless, theoretical 
estimates for the vacuum density are many orders of magnitude 
larger than its observed value. In addition, Da and D ln are of the 
same order of magnitude only at present, which marks our epoch 
as a special time in the evolution of the Universe (the “coinci¬ 
dence problem”). This lack of a clear theoretical understanding 
has motivated the development of a wide variety of alternative 
models. Those models which are close to ACDM are in broad 
agreement with current constraints on the background cosmol¬ 
ogy, but the perturbations may still evolve differently, and hence 
it is important to test their predictions against CMB data. 

There are at least three difficulties we had to face within this 
paper. First, there appears to be a vast array of possibilities in 
the literature and no agreement yet in the scientific community 
on a comprehensive framework for discussing the landscape of 
models. A second complication is that robust constraints on DE 
come from a combination of different data sets working in con¬ 
cert. Hence we have to be careful in the choice of the data sets 

1 Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci¬ 
entific consortia funded by ESA member states and led by Principal 
Investigators from France and Italy, telescope reflectors provided 
through a collaboration between ESA and a scientific consortium led 
and funded by Denmark, and additional contributions from NASA 
(USA). 


so that we do not find apparent hints for non-standard models 
that are in fact due to systematic errors. A third area of concern 
is the fact that numerical codes available at present for DE and 
MG are not as well tested in these scenarios as for ACDM, es¬ 
pecially given the accuracy reached by the data. Furthermore, 
in some cases, we need to rely on stability routines that deserve 
further investigation to assure that they are not excluding more 
models than required. 

In order to navigate the range of modelling possibilities, we 
adopt the following three-part approach. 

1 . Background parameterizations. Here we consider only pa- 
rameterizations of background-level cosmological quanti¬ 
ties. Perturbations are always included, but their evolution 
depends only on the background. This set includes models 
involving expansions, parameterizations or principal compo¬ 
nent analyses of the equation of state w = p/p of a DE fluid 
with pressure p and energy density p. Early DE also belongs 
to this class. 

2. Perturbation parameterizations. Here the perturbations 
themselves are parameterized or modified explicitly, not only 
as a consequence of a change in background quantities. 
There are two main branches we consider: firstly, effective 
field theory for DE (EFT, e.g. Gubitosi et al. 2013), which 
has a clear theoretical motivation, since it includes all the¬ 
ories derived when accounting for all symmetry operators 
in the Lagrangian, written in unitary gauge, i.e. in terms of 
metric perturbations only. This is a very general classifica¬ 
tion that has the advantage of providing a broad overview 
of (at least) all universally coupled DE models. However, a 
clear disadvantage is that the number of free parameters is 
large and the constraints are consequently weak. Moreover, 
in currently available numerical codes one needs to rely on 
stability routines which are not fully tested and may discard 
more models than necessary. 

As a complementary approach, we include a more phe¬ 
nomenological class of models obtained by directly param¬ 
eterizing two independent functions of the gravitational po¬ 
tentials. This approach can in principle probe all degrees of 
freedom at the background and perturbation level (e.g. Kunz 
2012) and is easier to handle in numerical codes. While the 
connection to physical models is less obvious here than in 
EFT, this approach allows us to gain a more intuitive under¬ 
standing of the general constraining power of the data. 

3. Examples of particular models. Here we focus on a se¬ 
lection of theories that have already been discussed in the 
literature and are better understood theoretically; these can 
partly be considered as applications of previous cases for 
which the CMB constraints are more informative, because 
there is less freedom in any particular theory than in a more 
general one. 

The CMB is the cleanest probe of large scales, which are 
of particular interest for modifications to gravity. We will inves¬ 
tigate the constraints coming from Planck data in combination 
with other data sets, addressing strengths and potential weak¬ 
nesses of different analyses. Before describing in detail the mod¬ 
els and data sets that correspond to our requirements, in Sect. 2 
we first address the main question that motivates our paper, dis¬ 
cussing why CMB is relevant for DE. We then present the spe¬ 
cific model parameterizations in Sect. 3. The choice of data sets 
is discussed in detail in Sec. 4 before we present results in Sect. 5 
and discuss conclusions in Sect. 6. 
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2. Why is the CMB relevant for dark energy? 

The CMB anisotropies are largely generated at the last-scattering 
epoch, and hence can be used to pin down the theory at early 
times. In fact many forecasts of future DE or MG experiments 
are for new data plus constraints from Planck. However, there 
are also several effects that DE and MG models can have on the 
CMB, some of which are to: 

1. change the expansion history and hence distance to the last 
scattering surface, with a shift in the peaks, sometimes re¬ 
ferred to as a geometrical projection effect (Hu & White 
1996); 

2. cause the decay of gravitational potentials at late times, af¬ 
fecting the low-multipole CMB anisotropies through the in¬ 
tegrated Sachs-Wolfe (or ISW) effect (Sachs & Wolfe 1967; 
Kofman & Starobinskii 1985); 

3. enhance the cross-correlation between the CMB and large- 
scale structure, through the ISW effect (Giannantonio et al. 
2008); 

4. change the lensing potential, through additional DE pertur¬ 
bations or modifications of GR (Acquaviva & Baccigalupi 
2006; Carbone et al. 2013); 

5. change the growth of structure (Peebles 1984; 
Barrow & Saich 1993) leading to a mismatch between the 
CMB-inferred amplitude of the fluctuations A s and late-time 
measurements of erg (Kunz et al. 2004; Baldi & Pettorino 
2011 ); 

6. impact small scales, modifying the damping tail in C T ( T , giv¬ 
ing a measurement of the abundance of DE at different red- 
shifts (Calabrese et al. 2011; Reichardt et al. 2012); 

7. affect the ratio between odd and even peaks if modifica¬ 
tions of gravity treat baryons and cold dark matter differently 
(Amendola et al. 2012b); 

8. modify the lensing B-mode contribution, through changes in 
the lensing potential (Amendola et al. 2014); 

9. modify the primordial B-mode amplitude and scale depen¬ 
dence, by changing the sound speed of gravitational waves 
(Antolini et al. 2013; Amendola et al. 2014; Raveri et al. 
2014). 

In this paper we restrict our analysis to scalar perturbations. 
The dominant effects on the temperature power spectrum are 
due to lensing and the ISW effect, as can be seen in Fig. 1, 
which shows typical power spectra of temperature anisotropies 
and lensing potential for modified gravity models. Different 
curves correspond to different choices of the p and // functions, 
which change the relation between the metric potentials and the 
sources, as well as introducing a gravitational slip; we will de¬ 
fine these functions in Sect. 3.2.2, Eq. (4) and Eq. (6), respec¬ 
tively. Spectra are obtained using a scale-independent evolution 
for both p and The two parameters in the figure then determine 
the change in amplitude of p and i] with respect to the ACDM 
case, in which E\\ = E 22 = 0 and p = rj = 1. 

3. Models and parameterizations 

We now provide an overview of the models addressed in this 
paper. Details on the specific parameterizations will be discussed 
in Sect. 5, where we also present the results for each specific 
method. 

We start by noticing that one can generally follow two dif¬ 
ferent approaches: (1) given a theoretical set up, one can specify 
the action (or Lagrangian) of the theory and derive background 



10 1 10 2 10 3 
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Fig. 1 . Typical effects of modified gravity on theoretical CMB 
temperature (top panel) and lensing potential (bottom panel) 
power spectra. An increase (or decrease) of E 22 with respect 
to zero introduces a gravitational slip, higher at present, when 
Q dc is higher (see Eq. (4) and Eq. (6)); this in turns changes the 
Weyl potential and leads to a higher (or lower) lensing poten¬ 
tial. On the other hand, whenever E\\ and E 22 are different from 
zero (quite independently of their sign) p and rj change in time: 
as the dynamics in the gravitational potential is increased, this 
leads to an enhancement in the ISW effect. Note also that even 
when the temperature spectrum is very close to ACDM (as for 
E 11 = E 22 = 0.5) the lensing potential is still different with re¬ 
spect to ACDM, shown in black. 


and perturbation equations in that framework; or (2) more phe¬ 
nomenologically, one can construct functions that map closely 
onto cosmological observables, probing the geometry of space- 
time and the growth of perturbations. Assuming spatial flatness 
for simplicity, the geometry is given by the expansion rate H and 
perturbations to the metric. If we consider only scalar-type com¬ 
ponents the metric perturbations can be written in terms of the 
gravitational potentials <E> and 'T (or equivalently by any two in¬ 
dependent combinations of these potentials). Cosmological ob¬ 
servations thus constrain one “background" function of time 
H(a) and two “perturbation" functions of scale and time a) 
and TBk, a) (e.g., Kunz 2012). These functions fix the metric, 
and thus the Einstein tensor G MV . Einstein’s equations link this 
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tensor to the energy-momentum tensor 7) JV , which in turn can be 
related to DE or MG properties. 

Throughout this paper we will adopt the metric given by the 
line element 

ds 2 = a 2 [-(1 + 2TOc/t 2 + (1 - 2€>)c/.y 2 J . (1) 

The gauge invariant potentials and 'F are related to 
the Bardeen (1980) potentials <I> A and <t> H and to the 
Kodama & Sasaki (1984) potentials v F K s and <E>ks in the follow¬ 
ing way: T = (Pa = 'F K s and <1> = -<1 >h = -®ks- Throughout 
the paper we use a metric signature (-, + + +) and follow the 
notation of Ma & Bertschinger (1995); the speed of light is set 
to c = 1, except where explicitly stated otherwise. 

We define the equation of state p(a) = w(a)p(a ), where p 
and p are the average pressure and energy density. The sound 
speed c s is defined in the fluid rest frame in terms of pres¬ 
sure and density perturbations as 6p{k, a) = c 2 (k, a)5p(k, a). The 
anisotropic stress cr(k,a) (equivalent to 7 Tt in the notation of 
Kodama & Sasaki 1984) is the scalar part of the off-diagonal 
space-space stress energy tensor perturbation. The set of func¬ 
tions {H, ®, V T) describing the metric is formally equivalent to 
the set of functions {w, , cr ) (Ballesteros et al. 2012). 

Specific theories typically cover only subsets of this function 
space and thus make specific predictions for their form. In the 
following sections we will discuss the particular theories that we 
consider in this paper. 

3.1. Background parameterizations 

The first main ‘category’ of theories we describe includes pa¬ 
rameterizations of background quantities. Even when we are 
only interested in constraints on background parameters, we 
are implicitly assuming a prescription for Dark Energy fluctua¬ 
tions. The conventional approach, that we adopt also here, is to 
choose a minimally-coupled scalar field model (Wetterich 1988; 
Ratra & Peebles 1988), also known as quintessence, which cor¬ 
responds to the choice of a rest-frame sound speed = 1 (i.e., 
equal to the speed of light) and cr - 0 (no scalar anisotropic 
stress). In this case the relativistic sound speed suppresses the 
dark energy perturbations on sub-horizon scales, preventing it 
from contributing significantly to clustering. 

Background parameterizations discussed in this paper in¬ 
clude: 

- (wo, w fl ) Taylor expansion at first order (and potentially 
higher orders); 

- Principal Component Analysis of w{a) (Huterer & Starkman 
2003), that allows to estimate constraints on w in indepen¬ 
dent redshift bins; 

- general parameterization of any minimally coupled scalar 
field in terms of three parameters e s , £ s , e°o. This is a 
novel way to describe minimally coupled scalar field mod¬ 
els without explicitly specifying the form of the potential 
(Huang etal. 2011); 

- Dark Energy density as a function of z (including parameter¬ 
izations such as early Dark Energy). 

The specific implementation for each of them is discussed 
in Sect. 5.1 together with corresponding results. We will con¬ 
clude the background investigation by describing, in Sect. 5.1.6, 
a compressed Gaussian likelihood that captures most of the 
constraining power of the Planck data applied to smooth Dark 
Energy or curved models (following Mukherjee et al. 2008). The 
compressed likelihood is useful for example to include more 


easily the Planck CMB data in Fisher-forecasts for future large- 
scale structure surveys. 


3.2. Perturbation parameterizations 

Modified gravity models (in which gravity is modified with re¬ 
spect to GR) in general affect both the background and the per¬ 
turbation equations. In this subsection we go beyond background 
parameterizations and identify two different approaches to con¬ 
strain MG models, one more theoretically motivated and a sec¬ 
ond more phenomenological one. We will not embark on a full- 
scale survey of DE and MG models here, but refer the reader to 
e.g. Amendola et al. (2013) for more details. 


3.2.1. Modified gravity and effective field theory 


The first approach starts from a Lagrangian, derived from an ef¬ 
fective field theory (EFT) expansion (Cheung et al. 2008), dis¬ 
cussed in Creminelli et al. (2009) and Gubitosi et al. (2013) in 
the context of DE. Specifically, EFT describes the space of (uni¬ 
versally coupled) scalar field theories, with a Lagrangian written 
in unitary gauge that preserves isotropy and homogeneity at the 
background level, assumes the weak equivalence principle, and 
has only one extra dynamical field besides the matter fields con¬ 
ventionally considered in cosmology. The action reads: 


S = 




d 4 x^\ [1 + D(r)] R + A(r) - a 2 c(f)8g 


.00 


M 4 (t) 


'1 

2 

* 2 / 


■ (a 2 dg oo y - M\(T)2a 2 5g 00 8K £ 

Mi{f) , \2 M 2 (t) o 2 M 2 (t) m 

(<«£) - -^-SK^SK; + -——-8g 00 8R (3> 

+ m\(T) (g^ + n IJ n v ) d M ( a 2 g m ) d v (rrg 00 ) j 
S m \xi> g/jv I - (2) 


Here R is the Ricci scalar, dR ( V> is its spatial perturbation, K 1 '. 
is the extrinsic curvature, and niQ is the bare (reduced) Planck 
mass. The matter part of the action, S m , includes all fluid 
components except dark energy, i.e., baryons, cold dark mat¬ 
ter, radiation, and neutrinos. The action in Eq. (2) depends 
on nine time-dependent functions (Bloomfield et al. 2013), here 
{Q.,c,K,M 2 ,M 4 ,M 2 ,M 4 ,M 2 ,m^}, whose choice specifies the 
theory. In this way, EFT provides a direct link to any scalar field 
theory. A particular subset of EFT theories are the Horndeski 
(1974) models, which include (almost) all stable scalar-tensor 
theories, universally coupled to gravity, with second-order equa¬ 
tions of motion in the fields and depend on five functions of 
time (Gleyzes et al. 2013; Bellini & Sawicki 2014; Piazza et al. 
2014). 

Although the EFT approach has the advantage of being very 
versatile, in practice it is necessary to choose suitable parameter¬ 
izations for the free functions listed above, in order to compare 
the action with the data. We will describe our specific choices, 
together with results for each of them, in Sect. 5.2. 


3.2.2. MG and phenomenological parameterizations 

The second approach adopted in this paper to test MG is more 
phenomenological and starts from the consideration that cosmo¬ 
logical observations probe quantities related to the metric pertur- 
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bations, in addition to the expansion rate. Given the line element 
of Eq. (1), the metric perturbations are determined by the two 
potentials <h and V P, so that we can model all observationally 
relevant degrees of freedom by parameterizing these two poten¬ 
tials (or, equivalently, two independent combinations of them) as 
functions of time and scale. Since a non-vanishing anisotropic 
stress (proportional to < 1 > - TO is a generic signature of modi¬ 
fications of GR (Mukhanov et al. 1992; Saltas et al. 2014), the 
parameterized potentials will correspond to predictions of MG 
models. 

Various parameterizations have been considered in the liter¬ 
ature. Some of the more popular (in longitudinal gauge) are: 

1. Q{a,k), which modifies the relativistic Poisson equation 
through extra DE clustering according to 

- k 2 d> = 4m Ga 2 Q(a, k)pA, (3) 

where A is the comoving density perturbation; 

2. iu(a,k) (sometimes also called Y(a, k)), which modifies the 
equivalent equation for V P rather than <i>: 

- k lx Y = 4nGa 2 p(a, k) P A\ (4) 

3. T.(a,k), which modifies lensing (with the lensing/Weyl po¬ 
tential being <t> + TO, such that 

- k 2 (d> +'¥) = 8nGa 2 T.(a, k)pA\ (5) 

4. 77 ( 72 , k), which reflects the presence of a non-zero anisotropic 
stress, the difference between d> and 'F being equivalently 
written as a deviation of the ratio 2 

77 (a, k) = <D/'P. ( 6 ) 

In the equations above, pA = p m A m + p, A, so that the parameters 
Q, p, or X quantify the deviation of the gravitational potentials 
from the value expected in GR due to perturbations of matter 
and relativistic particles. At low redshifts, where most DE mo¬ 
dels become relevant, we can neglect the relativistic contribu¬ 
tion. The same is true for 77 , where we can neglect the contribu¬ 
tion of relativistic particles to the anisotropic stress at late times. 

The four functions above are certainly not independent. It is 
sufficient to choose two independent functions of time and scale 
to describe all modifications with respect to General Relativity 
(e.g. Zhang et al. 2007; Amendola et al. 2008b). Popular choices 
include: (p, 77 ), which have a simple functional form for many 
theories; (p, X), which is more closely related to what we actu¬ 
ally observe, given that CMB lensing, weak galaxy lensing and 
the IS W effect measure a projection or derivative of the Weyl po¬ 
tential <t> + TO Furthermore, redshift space distortions constrain 
the velocity field, which is linked to 'P through the Euler equa¬ 
tion of motion. 

All four quantities, Q, p, X, and 77 , are free functions of time 
and scale. Their parameterization in terms of the scale factor a 
and momentum k will be specified in Section 5.2.2, together 
with results obtained by confronting this class of models with 
data. 

3.3. Examples of particular models 

The last approach is to consider particular models. Even 
though these are in principle included in the case described in 
Sect. 3.2.1, it is nevertheless still useful to highlight some well 
known examples of specific interest, which we list below. 

2 This parameter is called y in the code MGCAMB, but since y is also 
often used for the growth index, we prefer to use the symbol rj. 


- Minimally-coupled models beyond simple quintessence. 
Specifically, we consider “k-essence” models, which are 
defined by an arbitrary sound speed c 2 in addition to a 
free equation of state parameter w (Armendariz-Picon et al. 
2000 ). 

- An example of a generalized scalar field model 
(Deffayet et al. 2010) and of Lorentz-violating massive 
gravity (Dubovsky 2004; Rubakov & Tinyakov 2008), both 
in the ‘equation of state’ formalism of Battye & Pearson 
( 2012 ). 

- Universal “fifth forces.” We will show results for f(R) the¬ 
ories (Wetterich 1995a; Capozziello 2002; Amendola et al. 
2007; De Felice & Tsujikawa 2010), which form a subset of 
all models contained in the EFT approach. 

- Non-universal fifth forces. We will illustrate results for cou¬ 
pled DE models (Amendola 2000), in which dark matter par¬ 
ticles feel a force mediated by the DE scalar field. 

All these particular models are based on specific ac¬ 
tions, ensuring full internal consistency. The reviews by 
Amendola et al. (2013), Clifton et al. (2012), Joyce et al. (2014) 
and Huterer et al. (2015) contain detailed descriptions of a large 
number of models discussed in the literature. 

4. Data 

We now discuss the data sets we use, both from Planck and in 
combination with other experiments. As mentioned earlier, if 
we combine many different data sets (not all of which will be 
equally reliable) and take them all at face value, we risk attribut¬ 
ing systematic problems between data sets to genuine physical 
effects in DE or MG models. On the other hand, we need to avoid 
bias in confirming ACDM, and remain open to the possibility 
that some tensions may be providing hints that point towards DE 
or MG models. While discussing results in Sect. 5, we will try to 
assess the impact of additional data sets, separating them from 
the Planck baseline choice, keeping in mind caveats that might 
appear when considering some of them. For a more detailed dis¬ 
cussion of the data sets we refer to Planck Collaboration XIII 
(2015). 

4. 1. Planck data sets 

4.1.1. Planck low-f data 

The 2013 papers used WMAP polarization measurements 
(Bennett et al. 2013) at multipoles ( < 23 to constrain the op¬ 
tical depth parameter r. The corresponding likelihood was de¬ 
noted “WP” in the 2013 papers. 

For the present release, we use in its place a Planck 
polarization likelihood that is built through low-resolution 
maps of Stokes Q and U polarization measured by LFI at 
70 GHz (excluding data from Surveys 2 and 4), foreground- 
cleaned with the LFI 30 GHz and HFI 353 GHz maps, used 
as polarized synchrotron and dust templates, respectively (see 
Planck Collaboration XI (2015)). 

The foreground-cleaned LFI 70 GHz polarization maps 
are processed, together with the temperature map from the 
Commander component separation algorithm over 94% of the 
sky (see Planck Collaboration IX 2015, for further details), us¬ 
ing the low-f Planck temperature-polarization likelihood. This 
likelihood is pixel-based, extends up to multipoles £ = 29 and 
masks the polarization maps with a specific polarization mask, 
which uses 46% of the sky. Use of this likelihood is denoted as 
“lowP” hereafter. 
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The Planck lowP likelihood, when combined with the high- 
C Planck temperature one, provides a best fit value for the 
optical depth r = 0.078 + 0.019, which is about 1 cr lower 
than the value inferred from the WP polarization likelihood, 
i.e., r = 0.089 + 0.013, in the Planck 2013 papers (see also 
Planck Collaboration XIII 2015). However, we find that the LFI 
70 GHz and WMAP polarization maps are extremely consistent 
when both are cleaned with the HFI 353 GHz polarized dust 
template, as discussed in more detail in Planck Collaboration XI 
(2015). 


4.1.2. Planck high-f’data 

Following Planck Collaboration XV (2014), the high-f part of 
the likelihood (30 < ( < 2500) uses a Gaussian approximation, 

-log uc\cm = ^(C-C(60) t -C-' -(C-Cm + const. , (7) 


with C the data vector, C(6) the model with parameters 8 and C 
the covariance matrix. The data vector consists of the tempera¬ 
ture power spectra of the best CMB frequencies of the HFI in¬ 
strument. Specifically, as discussed in Planck Collaboration XI 
(2015), we use 100 GHz, 143 GHz and 217 GHz half-mission 
cross-spectra, measured on the cleanest part of the sky, avoid¬ 
ing the Galactic plane, as well as the brightest point sources 
and regions where the CO emission is the strongest. The point 
source masks are specific to each frequency. We retain, 66% 
of the sky for the 100 GHz map, 57 % for 143 GHz, and 47 % 
for 217 GHz. All the spectra are corrected for beam and pixel 
window functions. Not all cross-spectra and multipoles are in¬ 
cluded in the data vector; specifically, the TT 100 x 143 and 
100 x 217 cross-spectra, which do not bring much extra infor¬ 
mation, are discarded. Similarly, we only use multipoles in the 
range 30 < l < 1200 for 100 x 100 and 30 < £ < 2000 for 
143 x 143, discarding modes where the S/N is too low. We do 
not co-add the different cross-frequency spectra, since, even af¬ 
ter masking the highest dust-contaminated regions, each cross¬ 
frequency spectrum has a different, frequency-dependent resid¬ 
ual foreground contamination, which we deal with in the model 
part of the likelihood function. 

The model, C(8) can be rewritten as 


C,,M 


C cmb + C%{8) 

-yjA^Ay 


( 8 ) 


where C cmb is the set of CMB Qs, which is independent of 
frequency, C^ V (0) is the foreground model contribution to the 
cross-frequency spectrum /r x v, and A f , the calibration factor 
for the /r x /r spectrum. We retain the following contributions in 
our foreground modelling: dust; clustered cosmic infrared back¬ 
ground (CIB); thermal Sunyaev-Zeldovich (tSZ) effect; kinetic 
Sunyaev-Zeldovich (kSZ) effect; tSZ-CIB cross-correlations; 
and point sources. The dust, CIB and point source contributions 
are the dominant contamination. Specifically, dust is the domi¬ 
nant foreground at i < 500, while the diffuse point source term 
(and CIB for the 217 x 217) dominates the small scales. All our 
foreground models are based upon smooth Q templates with 
free amplitudes. All templates but the dust are based on analyti¬ 
cal models, as described in Planck Collaboration XI (2015). The 
dust is based on a mask difference of the 545 GHz map and is 
well described by a power law of index n = -2.63, with a wide 
bump around l = 200. A prior for the dust amplitude is com¬ 
puted from the cross-spectra with the 545 GHz map. We refer 


the reader to Planck Collaboration XI (2015) for a complete de¬ 
scription of the foreground model. The overall calibration for the 
100 x 100 and 217x217 power spectra free to vary within a prior 
measured on a small fraction of the sky near the Galactic pole. 

The covariance matrix C accounts for the correlations 
due to the mask and is computed following the equations in 
Planck Collaboration XV (2014). The fiducial model used to 
compute the covariance is based on a joint fit of ACDM and nui¬ 
sance parameters. The covariance includes the non-Gaussianity 
of the noise, but assumes Gaussian statistics for the dust. The 
non-whiteness of the noise is estimated from the difference be¬ 
tween the cross- and auto-half mission spectra and accounted 
for in an approximate manner in the covariance. Different Monte 
Carlo based corrections are applied to the covariance matrix cal¬ 
culation to account for inaccuracies in the analytic formulae at 
large scales (l < 50) and when dealing with the point source 
mask. Beam-shape uncertainties are folded into the covariance 
matrix. A complete description of the computation and its vali¬ 
dation is discussed in Planck Collaboration XI (2015). 

The TT unbinned covariance matrix is of size about 8000 x 
8000. When adding the polarization, the matrix has size 23000 x 
23000, which translates into a significant memory requirement 
and slows the likelihood computation considerably. We thus bin 
the data and covariance matrix, using a variable bin-size scheme, 
to reduce the data vector dimension by about a factor of ten. We 
checked that for the ACDM model, including single parameter 
classical extensions, the cosmological and nuisance parameter 
fits are identical with or without binning. 


4.1.3. Planck CMB lensing 

Gravitational lensing by large-scale structure introduces depen¬ 
dencies in CMB observables on the late-time geometry and clus¬ 
tering, which otherwise would be degenerate in the primary 
anisotropies (Hu 2002; Lewis & Challinor 2006). This provides 
some sensitivity to dark energy and late-time modifications of 
gravity from the CMB alone. The source plane for CMB lensing 
is the last-scattering surface, so the peak sensitivity is to lenses 
at z ~ 2 (i.e., half-way to the last-scattering surface) with typi¬ 
cal sizes of order 10 2 Mpc. Although this peak lensing redshift is 
rather high for constraining simple late-time dark energy mod¬ 
els, CMB lensing deflections at angular multipoles l ^ 60 have 
sources extending to low enough redshift that DE becomes dy¬ 
namically important (e.g., Pan et al. 2014). 

The main observable effects of CMB lensing are a smooth¬ 
ing of the acoustic peaks and troughs in the temperature 
and polarization power spectra, the generation of significant 
non-Gaussianity in the form of a non-zero connected 4-point 
function, and the conversion of E-mode to fi-mode polar¬ 
ization. The smoothing effect on the power spectra is in¬ 
cluded routinely in all results in this paper. We addition¬ 
ally include measurements of the power spectrum of the 
CMB lensing potential 4k which are extracted from the Planck 
temperature and polarization 4-point functions, as presented 
in Planck Collaboration XV (2015) and discussed further below. 
Lensing also produces 3-point non-Gaussianity, which peaks 
in squeezed configurations, due to the correlation between the 
lensing potential and the ISW effect in the large-angle tem¬ 
perature anisotropies. This effect has been measured at around 
3 cr with the full-mission Planck data (Planck Collaboration XV 
2015; Planck Collaboration XXI 2015). Although in principle 
this is a further probe of DE (Verde & Spergel 2002) and 
MG (Acquaviva et al. 2004), we do not include these T-cj) corre- 
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lations in this paper as the likelihood was not readily available. 
We plan however to test this effect in future work. 

The construction of the CMB lensing likelihood we use 
in this paper is described fully in Planck Collaboration XV 
(2015); see also Planck Collaboration XIII (2015). It is a sim¬ 
ple Gaussian approximation in the estimated bandpow- 

ers, covering the multipole range 40 < l < 400. The are 
estimated from the full-mission temperature and polariza¬ 
tion 4-point functions, using the SMICA component-separated 
maps (Planck Collaboration IX 2015) over approximately 70 % 
of the sky. A large number of tests of internal consistency 
of the estimated to different data cuts (e.g., whether 
polarization is included, or whether individual frequency 
bands are used in place of the SMICA maps) are reported 
in Planck Collaboration XV (2015). All such tests are passed 
for the conservative multipole range 40 < ( < 400 that we 
adopt in this paper. For multipoles l > 400, there is marginal 
evidence of systematic effects in reconstructions of the lens¬ 
ing deflections from temperature anisotropies alone, based on 
curl-mode tests. Reconstructing the lensing deflections on large 
angular scales is very challenging because of the large “mean- 
field” due to survey anisotropies, which must be carefully sub¬ 
tracted with simulations. We conservatively adopt a minimum 
multipole of i - 40 here, although the results of the null tests 
considered in Planck Collaboration XV (2015) suggest that this 
could be extended down to l = 8. For Planck, the multipole 
range 40 < ( < 400 captures the majority of the S/N on Cp‘ 
for ACDM models, although this restriction may be more lossy 
in extended models. The Planck 2014 lensing measurements are 
the most significant to date (the amplitude of C'y’ is measured at 
greater than 40 cr), and we therefore choose not to include lens¬ 
ing results from other CMB experiments in this paper. 


In this paper we will only show results that include CMB 
high-£ polarization data where we find that it has a signif¬ 
icant impact. DE and MG can in principle also affect the 
B-mode power spectrum through lensing of B-modes (if the 
lensing Weyl potential is modified) or by changing the po¬ 
sition and amplitude of the primordial peak (Antolini et al. 
2013; Pettorino & Amendola 2014), including modifications of 
the sound speed of gravitational waves (Amendola et al. 2014; 
Raveri et al. 2014). Due to the unavailability of the likelihood, 
results from B-mode polarization are left to future work. 

4.2. Background data combination 

We identify a first basic combination of data sets that we mostly 
rely on, for which we have a high confidence that systematics are 
under control. Throughout this paper, we indicate for simplicity 
with “BSH” the combination BAO + SN-Ia + Hq, which we now 
discuss in detail. 

4.2.1. Baryon acoustic oscillations 

Baryon acoustic oscillations (BAO) are the imprint of oscilla¬ 
tions in the baryon-photon plasma on the matter power spec¬ 
trum and can be used as a standard ruler, calibrated to the CMB- 
determined sound horizon at the end of the drag epoch. Since 
the acoustic scale is so large, BAO are largely unaffected by 
nonlinear evolution. As in the cosmological parameter paper, 
Planck Collaboration XIII (2015), BAO is considered as the pri¬ 
mary data set to break parameter degeneracies from CMB mea¬ 
surements and offers constraints on the background evolution of 
MG and DE models. The BAO data can be used to measure both 
the angular diameter distance Z>a(z), and the expansion rate of 
the Universe H(z) either separately or through the combination 


4.1.4. Planck CMB polarization 

The TE and EE likelihood follows the same principle as the TT 
likelihood described in Sect. 4.1.2. The data vector is extended 
to contain the TE and EE cross-half-mission power spectra 
of the same 100 GHz, 143 GHz, and 217 GHz frequency maps. 
Following Planck Collaboration Int. XXX (2014), we mask the 
regions where the dust intensity is important, and retain 70 %, 
50%, and 41 % of the sky for our three frequencies. We ig¬ 
nore any other polarized galactic emission and in particular 
synchrotron, which has been shown to be negligible, even at 
100 GHz. We use all of the cross-frequency spectra, using the 
multipole range 30 < l < 1000 for the 100 GHz cross-spectra 
and 500 < l < 2000 for the 217 GHz cross-spectra. Only the 
143 x 143 spectrum covers the full 30 < £ < 2000 range. 
We use the same beams as for the TT spectra and do not cor¬ 
rect for leakage due to beam mismatch. A complete descrip¬ 
tion of the beam mismatch effects and correction is described 
in Planck Collaboration XI (2015). 

The model is similar to the TT one. We retain a single fore¬ 
ground component accounting for the polarized emission of the 
dust. Following Planck Collaboration Int. XXX (2014), the dust 
C( template is a power law with index n = -2.4. A prior for 
the dust amplitude is measured in the cross-correlation with the 
353 GHz maps. The calibration parameters are fixed to unity. 

The covariance matrix is extended to polarization, as de¬ 
scribed in Planck Collaboration XI (2015), using the correlation 
between the TT, TE, and EE spectra. It is computed similarly 
to the TT covariance matrix, as described in Sect. 4.1.2. 


Dy(z) 


(1 +z) 2 DX(z) 


cz 

HU) 


1/3 


(9) 


As in Planck Collaboration XIII (2015) we use: the SDSS 
Main Galaxy Sample at z e ff = 0.15 (Ross et al. 2014); the 
Baryon Oscillation Spectroscopic Survey (BOSS) “LOWZ” 
sample at z e ft = 0.32 (Anderson et al. 2014); the BOSS CMASS 
(i.e. “constant mass” sample) at z c ff = 0.57 of Anderson et al. 
(2014); and the six-degree-Field Galaxy survey (6dFGS) at 
z e ff = 0.106 (Beutler et al. 2011). The first two measurements 
are based on peculiar velocity field reconstructions to sharpen 
the BAO feature and reduce the errors on the quantity Dy/r s ; the 
analysis in Anderson et al. (2014) provides constraints on both 
E>A(z e ff) and H(z e «). In all cases considered here the BAO obser¬ 
vations are modelled as distance ratios, and therefore provide no 
direct measurement of Hq. However, they provide a link between 
the expansion rate at low redshift and the constraints placed by 
Planck at z ~ 1100. 


4.2.2. Supernovae 

Type-la supernovae (SNe) are among the most important probes 
of expansion and historically led to the general acceptance that 
a DE component is needed (Riess et al. 1998; Perlmutter et al. 
1999). Supernovae are considered as “standardizable candles” 
and so provide a measurement of the luminosity distance as a 
function of redshift. However, the absolute luminosity of SNe 
is considered uncertain and is marginalized out, which also re¬ 
moves any constraints on Hq. 
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Consistently with Planck Collaboration XIII (2015), we use 
here the analysis by Betoule et al. (2013) of the “Joint Light- 
curve Analysis” (JLA) sample. JLA is constructed from the 
SNLS and SDSS SNe data, together with several samples of 
low redshift SNe. Cosmological constraints from the JLA sam¬ 
ple 1 are discussed by Betoule et al. (2014), and as mentioned in 
Planck Collaboration XIII (2015) the constraints are consistent 
with the 2013 and 2104 Planck values for standard ACDM. 

4.2.3. The Hubble constant 

The CMB measures mostly physics at the epoch of recombina¬ 
tion, and so provides only weak direct constraints about low- 
redshift quantities through the integrated Sachs-Wolfe effect and 
CMB lensing. The CMB-inferred constraints on the local ex¬ 
pansion rate Hq are model dependent, and this makes the com¬ 
parison to direct measurements interesting, since any mismatch 
could be evidence of new physics. 

Here, we rely on the re-analysis of the Riess et al. (2011) 
(hereafter Rll) Cepheid data made by Efstathiou (2014) (here¬ 
after E14). By using a revised geometric maser distance to NGC 
4258 from Humphreys et al. (2013), E14 obtains the following 
value for the Hubble constant: 

H 0 = (70.6 ± 3.3)kins -1 Mpc' 1 , (10) 

which is within 1 cr of the Planck TT+lowP estimate. In this pa¬ 
per we use Eq. (10) as a conservative Ho prior. We note that 
the 2015 Planck TT+lowP value is perfectly consistent with 
the 2013 Planck value (Planck Collaboration XVI 2014) and 
so the tension with the Rll Ho determination is still present 
at about 2.4 cr. We refer to the cosmological parameter paper 
Planck Collaboration XIII (2015) for a more comprehensive dis¬ 
cussion of the different values of Hq present in the literature. 

4.3. Perturbation data sets 


In linear theory, anisotropic clustering along the line of 
sight and in the transverse directions measures the combination 
/(z)crg(z), where the growth rate is defined by 


m = 


dlncrg 
din a 


( 11 ) 


where cr 8 is calculated including all matter and neutrino den¬ 
sity perturbations. Anisotropic clustering also contains ge¬ 
ometric information from the Alcock-Paczynski (AP) ef¬ 
fect (Alcock & Paczynski 1979), which is sensitive to 


F ap (z) = (1+z)D a (z)H(z). (12) 


In addition, fits which constrain RSD frequently also mea¬ 
sure the BAO scale, Dy(z)/r s , where r s is the comoving sound 
horizon at the drag epoch, and D v is given in Eq. (9). As 
in Planck Collaboration XIII (2015) we consider only analyses 
which solve simultaneously for the acoustic scale, T/p and /erg. 

The Baryon Oscillation Spectroscopic Survey (BOSS) col¬ 
laboration have measured the power spectrum of their CMASS 
galaxy sample (Beutler et al. 2014) in the range k = 0.01- 
0.20 hMpW 1 . Samushia et al. (2014) have estimated the mul¬ 
tipole moments of the redshift-space correlation function of 
CMASS galaxies on scales > 25 h 'Mpc. Both papers provide 
tight constraints on the quantity /cr 8 , and the constraints are 
consistent. The Samushia et al. (2014) result was shown to be¬ 
have marginally better in terms of small-scale bias compared 
to mock simulations, so we choose to adopt this as our base¬ 
line result. Note that when we use the data of Samushia et al. 
(2014), we exclude the measurement of the BAO scale, Dy/r s , 
from Anderson et al. (2013), to avoid double counting. 

The Samushia et al. (2014) results are expressed as a 3 x 3 
covariance matrix for the three parameters Dv/r s , Tap and 
/erg, evaluated at an effective redshift of Zdr = 0.57. Since 
Samushia et al. (2014) do not apply a density field reconstruc¬ 
tion in their analysis, the BAO constraints are slightly weaker 
than, though consistent with, those of Anderson et al. (2014). 


The additional freedom present in MG models can be calibrated 
using external data that test perturbations in particular. In the 
following we describe other available data sets that we included 
in the grid of runs for this paper. 

4.3.1. Redshift space distortions 

Observations of the anisotropic clustering of galaxies in red¬ 
shift space permit the measurement of their peculiar velocities, 
which are related to the Newtonian potential 'F via the Euler 
equation. This, in turn, allows us to break a degeneracy with 
gravitational lensing that is sensitive to the combination <t> + T. 
Galaxy redshift surveys now provide very precise constraints on 
redshift-space clustering. The difficulty in using these data is 
that much of the signal currently comes from scales where non¬ 
linear effects and galaxy bias are significant and must be accu¬ 
rately modelled (see, e.g., the discussions in Bianchi et al. 2012; 
Gil-Marfn et al. 2012). Moreover, adopting the wrong fiducial 
cosmological model to convert angles and redshifts into dis¬ 
tances can bias measurements of the rate-of-growth of structure 
(Reid et al. 2012; Howlett et al. 2014). Significant progress in 
the modelling has been achieved in the last few years, so we 
shall focus here on the most recent (and relatively conservative) 
studies. A compilation of earlier measurements can be found in 
the references above. 

3 A CosmoMC likelihood module for the JLA sample is available at 
http://supernovae.in2p3.fr/sdss_snls_jla/ReadMe.html. 


4.3.2. Galaxy weak lensing 

The distortion of the shapes of distant galaxies by large-scale 
structure along the line of sight (weak gravitational lensing or 
cosmic shear) is particularly important for constraining DE and 
MG, due to its dependence on the growth of fluctuations and the 
two scalar metric potentials. 

Currently the largest weak lensing (WL) survey is the 
Canada France Hawaii Telescope Lensing Survey (CFHTLenS), 
and we make use of two data sets from this survey: 

1. 2D CFHTLenS data (Kilbinger et al. 2013), whose shear 
correlation functions are estimated in the angular range 
0.9 to 296.5 aremin; 

2. the tomographic CFHTLenS blue galaxy sample 
(Heymans et al. 2013), whose data have an intrinsic 
alignment signal consistent with zero, eliminating the need 
to marginalize over any additional nuisance parameters, and 
where the shear correlation functions are estimated in six 
redshift bins, each with an angular range 1.7 < 8 < 37.9 
aremin. 

Since these data are not independent we do not combine them, 
but rather check the consistency of our results with each. The 
galaxy lensing convergence power spectrum, T’L(T), can be writ¬ 
ten in terms of the Weyl potential, Pq+v, by 

p K u (e) * 2 T?c J %, Or) gj Or) /w (e/ x ,x) , G3) 
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Fig. 2. Q m -<x 8 constraints for tomographic lensing from 
Heymans et al. (2013), using a very conservative angular cut, as 
described in the text (see Sect. 4.3.2). We show results using lin¬ 
ear theory, nonlinear corrections from Halofit (HL) versions 
1, 4, marginalization over baryonic AGN feedback (BF), and 
intrinsic alignment (IA) (the latter two using nonlinear correc¬ 
tions and Halofit 4). Coloured points indicate Ho values from 
WL+HL4. 


where we have made use of the Limber approximation in flat 
space, and x is the comoving distance. The lensing efficiency is 
given by 

X oo / 

4 x'rii (x') ~ , X , (14) 

X 

where «, (x) is the radial distribution of source galaxies in bin i. 
In the case of no anisotropic stress and no additional clustering 
from the DE, the convergence power spectrum can be written in 
the usual form 

P-fO = £° 8l(x J 8 ^ } P(tlx,x)dx • (15) 

However, in this paper we always use the full Weyl potential to 
compute the theoretical WL predictions. The convergence can 
also be written in terms of the correlation functions via 

fijW = ^ f d UP*j(£)J±m, (16) 

where the Bessel functions are J + = Jo and J = ,/ 4 . 

In this paper we need to be particularly careful about the 
contribution of nonlinear scales to since the behaviour of MG 
models in the nonlinear regime is not known very precisely. The 
standard approach is to correct the power spectrum on nonlinear 
scales using the Halofit fitting function. Since its inception, 
there have been several revisions to improve the agreement with 
/V-body simulations. We use the following convention to label 
the particular Halofit model: 

1. the original model of Smith et al. (2003); 

2. an update from higher resolution /V-body simulations, to in¬ 
clude the effect of massive neutrinos (Bird et al. 2012); 

3. an update to improve the accuracy on small scales 4 ; 

4 http://www.roe.ac.uk/~j ap/haloes/ 


4. an update from higher resolution /V-body simulations, in¬ 
cluding DE cosmologies with constant equation of state 

(Takahashi et al. 2012). 

Given this correction, one can scale the Weyl potential transfer 
functions by the ratio of the nonlinear to linear matter power 
spectrum 

lp™ nlin (k,z ) 

7W£, z) -» 7W*. z) yj • (17) 

Both (Kilbinger et al. 2013) and (Heymans et al. 2013) quote 
a “conservative” set of cuts to mitigate uncertainty over the non¬ 
linear modelling scheme. For the 2D analysis of Kilbinger et al. 
(2013) angular scales 8 < 17' are excluded for /' + , and 8 < 54' 
for £~. For the tomographic analysis of Heymans et al. (2013), 
angular scales 8 < 3' are excluded for /'*' for any bin combina¬ 
tion involving the two lowest redshift bins, and no cut is applied 
for the highest four redshift bins. For , angular scales 8 < 30' 
are excluded for any bin combination involving the four lowest 
redshift bins, and 8 < 16' for the highest two bins. 

These cuts, however, may be insufficient for our purposes, 
since we are interested in extensions to ACDM. We therefore 
choose a very conservative set of cuts to mitigate the total con¬ 
tribution from nonlinear scales. In order to select these cuts we 
choose the baseline Planck TT+lowP ACDM cosmology as de¬ 
scribed in Planck Collaboration XIII (2015), for which one can 
use Eq. (15). The cuts are then chosen by considering Ax 1 = 
\Xiin ~TnoniiJ °f t ^ le WL likelihood as a function of angular cut. 
In order for this to remain Ax 2 < 1 for each of the Halofit 
versions, we find it necessary to remove entirely from each 
data set, and exclude 8 < 17' for f 4 for both the 2D and tomo¬ 
graphic bins. We note that a similar approach to Kitching et al. 
(2014) could also be followed using 3D CFHTLenS data, where 
the choice of cut is more well defined in //-space, however the 
likelihood for this was not available at the time of this paper. 

On small scales, the effects of intrinsic alignments and bary¬ 
onic feedback can also become significant. In order to check 
the robustness of our cuts to these effects we adopt the same 
methodology of MacCrann et al. (2014). Using the same base¬ 
line model and choosing Halofit version 4, we scale the mat¬ 
ter power spectrum by an active galactic nuclei (AGN) com¬ 
ponent, derived from numerical simulations (van Daalen et al. 
2011), marginalizing over an amplitude ctagn- The AGN bary¬ 
onic feedback model has been shown by Harnois-Deraps et al. 
(2014) to provide the best fit to small-scale CFHTLens data. For 
intrinsic alignment we adopt the model of Bridle & King (2007), 
including the additional nonlinear alignment contributions to , 
and again marginalizing over an amplitude (Y\\. For more details 
on this procedure, we refer the reader to MacCrann et al. (2014). 

The robustness of our ultra-conservative cuts to nonlinear 
modelling, baryonic feedback and intrinsic alignment margina¬ 
lization, is illustrated in Fig. 2 for the tomographic data, with 
similar constraints obtained from 2D data. Assuming the same 
base ACDM cosmology, and applying priors of Q h /i 2 = 0.0223 + 
0.0009, n s = 0.96 + 0.02, and 40kms _1 Mpc _1 < Ho < 
lOOkms 1 Mpc 1 to avoid over-fitting the model, we find that 
the WL likelihood is insensitive to nonlinear physics. We there¬ 
fore choose to adopt the tomographic data with the ultra¬ 
conservative cuts as our baseline data set. 

4.4. Combining data sets 

We show for convenience in Table 1 the schematic summary 
of models. All models have been tested for the combina- 
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Table 1 . Table of models tested in this paper. We have tested all models for the combinations: Planck, Planck+ BSH, Planck+ WL, 
Planck+ BAO/RSD and fVancAr+WL+BAO/RSD. Throughout the text, unless otherwise specified, Planck refers to the baseline 
Planck TT+lowP combination. The effects of CMB lensing and Planck TT,TE,EE polarization have been tested on all runs above 
and are, in particular, used to constrain the amount of DE at early times. 


Model 

Section 

ACDM. 

Planck Collaboration XIII (2015) 

Background parameterizations: 

w. 

Wo, W a . 

w higher order expansion. 

1-parameter w(a ). 

w PCA. 

e s ,^,e m . 

Early DE. 

Planck Collaboration XIII (2015) 
Sect. 5. 1.1: Figs. 3, 4, 5 

Sect. 5.1.1 

Sect. 5.1.2: Fig. 6 

Sect. 5.1.3: Fig. 7 

Sect. 5.1.4: Figs. 8, 9 

Sect. 5.1.5: Figs. 10, 11 

Perturbation parameterizations: 

EFT exponential . 

EFT linear. 

p, 7] scale-independent: 

DE-related. 

time related. 

p, r] scale-dependent: . 

DE-related. 

time related. 

Sect. 5.2.1: Fig. 12 

Sect. 5.2.1: Fig. 13 

Sect. 5.2.2: Figs. 1, 14, 15, 16, 17 
Sect. 5.2.2: Figs. 14, 16 

Sect. 5.2.2: Fig. 18 

Sect. 5.2.2 

Other particular examples: 

DE sound speed and k-essence . . 
Equation of state approach: .... 
Lorentz-violating massive gravity 
Generalized scalar fields. 

m . 

Coupled DE. 

Sect. 5.3.1 

Sect. 5.3.2 

Sect. 5.3.2 

Sect. 5.3.3: Figs. 19, 20 

Sect. 5.3.4: Figs. 21, 22 


tions: Planck, Planck+ BSH, Planck+ WL, Planck+BAOfRSD 
and Planck+ WL+BAO/RSD. Throughout the text, unless other¬ 
wise specified, Planck refers to the baseline Planck TT+lowP 
combination. The effects of CMB lensing and Planck TT,TE,EE 
polarization have been tested on all runs above and are, in parti¬ 
cular, used to constrain the amount of DE at early times. For each 
of them we indicate the section in which the model is described 
and the corresponding figures. In addition, all combinations in 
the table have been tested with and without CMB lensing. The 
impact of Planck high-f’ polarization has been tested on all mo¬ 
dels for the combination Planck+BAO+SNe+H 0 . 

5 . Results 

We now proceed by illustrating in detail the models and parame- 
terizations described in Sect. 3, through presenting results for 
each of them. The structure of this section is as follows. We 
start in Sect. 5.1 with smooth dark energy models that are ef¬ 
fectively parameterized by the expansion history of the Universe 
alone. In Sect. 5.2 we study the constraints on the presence of 
non-negligible dark energy perturbations, both in the context 
of general modified gravity models described through effective 
field theories and with phenomenological parameterizations of 
the gravitational potentials and their combinations, as illustrated 
in Sect. 3.2.2. The last part. Sect. 5.3, illustrates results for a 
range of particular examples often considered in the literature. 

5.1. Background parameterizations 

In this section, we consider models where DE is a generic 
quintessence-like component with equation of state w = p/p. 


where p and p are the spatially averaged (background) DE pres¬ 
sure and density. Although it is important to include, as we do, 
DE perturbations, models in this section have a sound speed that 
is equal to the speed of light, which means that they are smooth 
on sub-horizon scales (see Sect. 3.1 for more details). We start 
with Taylor expansions and a principal component analysis of 
w in a fluid formalism, then consider actual quintessence mo¬ 
dels parameterized through their potentials and finally study the 
limits that can be put on the abundance of DE density at early 
times. At the end of the sub-section we provide the necessary in¬ 
formation to compress the Planck CMB power spectrum into a 
4-parameter Gaussian likelihood for applications where the full 
likelihood is too unwieldy. 

5.1.1. Taylor expansions of w and w 0 , w a parameterization 

If the dark energy is not a cosmological constant with w - — 1 
then there is no reason why w should remain constant. In order 
to test a time-varying equation of state, we expand w{a) in a 
Taylor series. The first order corresponds to the {wo,w fl } case, 
also discussed in Planck Collaboration XIII (2015): 

w(a) = w 0 + (1 - a)w a . (18) 

We use the parameterized post-Friedmann (PPF) model of 
Hu & Sawicki (2007) and Fang et al. (2008) to allow for val¬ 
ues w < - 1 (note that there is another PPF formalism dis¬ 
cussed in Baker et al. (2014a)). Marginalized posterior distribu¬ 
tions for wo, w a , //(i and cr 8 are shown in Fig. 3 and the cor¬ 
responding 2D contours can be found in Fig. 4 for w a vs wo 
and for erg vs Q m . Results from Planck TT+lowP+BSH data 
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- Planck TT+lowP - Planck TT+lowP+WL - Planck TT+lowP+WL+BAO/RSD 

Planck TT+lowP+BSH Planck TT+1 ow P.BA0/RSD 



Fig. 3. Parameterization {wo, w a } (see Sect. 5.1.1). Marginalized posterior distributions for wo, w a , Hq and erg for various data combi¬ 
nations. The tightest constraints come from the Planck TT+lowP+BSH combination, which indeed tests background observations, 
and is compatible with ACDM. 


1.20 


1.05 


b 

0.90 


0.75 


- 2-10 1 
W 0 

Fig. 4. Marginalized posterior distributions of the {wq, w a ) parameterization (see Sect. 5.1.1) for various data combinations. The best 
constraints come from the priority combination and are compatible with ACDM. The dashed lines indicate the point in parameter 
space (-1,0) corresponding to the ACDM model. CMB lensing and polarization do not significantly change the constraints. Here 
Planck indicates Planck TT+lowP. 




are shown in blue and corresponds to the combination we con¬ 
sider the most secure, which in this case also gives the strongest 
constraints. This is expected, since the BAO and SNe data in¬ 
cluded in the BSH combination provide the best constraints on 
the background expansion rate. Results for weak lensing (WL) 
and redshift space distortions (RSD) are also shown, both sep¬ 
arately and combined. The constraints from these probes are 
weaker, since we are considering a smooth dark energy model 
where the perturbations are suppressed on small scales. While 
the WL data appear to be in slight tension with ACDM, accord¬ 
ing to the green contours shown in Fig. 4, the difference in total 
X 2 between the best-fit in the {wo, w fl ) model and in ACDM for 
Planck TT+lowP+WL is Ay 2 = -5.6, which is not very signif¬ 
icant for 2 extra parameters (for normal errors a 2 cr deviation 
corresponds to a x 2 absolute difference of 6.2). The WL con¬ 
tributes a A;^ 2 of -2.0 and the Ay 2 MB = -3.3 (virtually the same 
as when using Planck TT+lowP alone, for which &Xcmb ~ _ 3.2, 
which seems to indicate that WL is not in tension with Planck 
TT+lowP within a (wo,w a ) cosmology). However, as also dis¬ 
cussed in Planck Collaboration XIII (2015), these data combina¬ 
tions prefer very high values of Ho, which is visible also in the 
third panel of Fig. 3. The combination Planck TT+lowP+BSH, 


on the other hand, is closer to ACDM, with a total x 2 differ¬ 
ence between (wo,w fl ) and ACDM of only -0.8. We also show 
in Fig. 5 the equation of state reconstructed as a function of red- 
shift from the linear expansion in the scale factor a for different 
combinations of data. 

One might wonder whether it is reasonable to stop at first 
order in w(a). We have therefore tested a generic expansion in 
powers of the scale factor up to order N: 


N 

w(a) = wp + J^(l - a)'wi. (19) 

i=i 


We find that all parameters are very stable when allowing higher 
order polynomials; the w, parameters are weakly constrained and 
going from N = 1 (the linear case) to N = 2 (quadratic case) to 
N = 3 (cubic expansion) does not improve the goodness of fit 
and stays compatible with ACDM, which indicates that a linear 
parameterization is sufficient. 
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5.1.2. 1 -parameter varying w 


A simple example of a varying w model that can be written in 
terms of one extra parameter only (instead of wq, w a ) was pro¬ 
posed in Gott & Slepian (201 1); Slepian et al. (2014), motivated 
in connection to a DE minimally-coupled scalar field, slowly 
rolling down a potential \m 2 (f) 2 , analogous to the one predicted 
in chaotic inflation (Linde 1983). More generally, one can fully 
characterize the background by expanding a varying equation of 
state w(z) = -1 + 6w(z) ~ -1 + Swq x H^/H 2 (z), where: 


2 ^ ~ ^m( 1 + zY + ^de 
H o 


a+zy 


nSwo/fide 


^m(l 'b zY "b ^de 


( 20 ) 


at first order in <5wo, which is then the only extra parameter. 
Marginalized posterior contours in the plane 1i-5wq are shown 
in Fig. 6. The tightest constraints come from the combination 
Planck TT+lowP+lensing+BSH that gives 6wq = -0.008 + 
0.068 at 68 % confidence level, which slightly improves con¬ 
straints found by Aubourg et al. (2014). 


5.1.3. Principal Component Analysis on w(z) 



Swq 


Fig. 6. Marginalized posterior contours in the h-6wo plane are 
shown for 1-parameter varying w models (see Sect. 5.1.2) 
for different data combinations. Here Planck indicates Planck 
TT+lowP. 


A complementary way to measure the evolution of the equation 
of state, which is better able to model rapid variations, proceeds 
by choosing w in N fixed bins in redshift and by performing a 
principal component analysis to uncorrelate the constraints. We 
consider A' = 4 different bins in z and assume that w has a con¬ 
stant value pi in each of them. We then smooth the transition 
from one bin to the other such that: 


w(z) = i + Aw tanh 


Z - Zi 


+ 1 forz < Zi, i e {1, 4), (21) 


with Aw = (pj - pi-i)/ 2, a smoothing scale s = 0.02, and a 
binning Zi = (0,0.2,0.4,0.6,1.8). We have tested also a larger 
number of bins (up to N = 18) and have found no improvement 
in the goodness of fit. 



Fig. 5. Reconstructed equation of state w(z) as a function of red- 
shift (see Sect. 5.1.1), when assuming a Taylor expansion of w(z) 
to first-order (N - 1 in Eq. 19), for different combinations of 
the data sets. The coloured areas show the regions which con¬ 
tain 95 % of the models. The central blue line is the median 
line for Planck TT+lowP+BSH. Here Planck indicates Planck 
TT+lowP. 


The constraints on the vector pSi = of values 

that w(z) can assume in each bin is difficult to interpret, due to 
the correlations between bins. To uncorrelate the bins, we per¬ 
form a principal component analysis (Huterer & Starkman 2003; 
Huterer & Cooray 2005; Said et al. 2013). We first run COSMOMC 
(Lewis & Bridle 2002) on the original binning values p,\ then 
extract the covariance matrix that refers to the parameters we 
want to constrain: 


C = (pp T ) - (p)(p T ), (22) 


where p is the vector of parameters /?, and p [ is its transpose. 
We calculate the Fisher matrix, F = C 1 , and diagonalize it, 
F = 0 [ \0, where A is diagonal and O is the orthogonal ma¬ 
trix whose columns are the eigenvectors of the Fisher matrix. 
We then define W = 0 T A 1/2 0 (e.g., Huterer & Cooray 2005) 
and normalize this such that its rows sum up to unity; this ma¬ 
trix can be used to find the new vector q = Wp of uncorrelated 
parameters that describe w(z). This choice of W has been shown 
to be convenient, since most of the weights (i.e., the rows of 
W ) are found to be positive and fairly well localized in redshift. 
In Fig. 7 (lower panel) we show the weights for each bin as a 
function of redshift. Because they overlap only partially, we can 
assume the binning to be the same as the original one and attach 
to each of them error bars corresponding to the mean and stan¬ 
dard deviations of the q values. The result is shown in Fig. 7, top 
panel. The equation of state is compatible with the ACDM value 
w — — 1. Note however that this plot contains more information 
than a Taylor expansion to first order. 


5.1.4. Parameterization for a weakly-coupled canonical 
scalar field. 

We continue our investigation of background parameterizations 
by considering a slowly rolling scalar field. In this case, as in 
inflation, we can avoid writing down an explicit potential V((f>) 
and instead parameterize w{a) at late times, in the presence of 
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0.0 0.5 1.0 1.5 2.0 



Fig. 7. PCA analysis constraints (described in Sect. 5.1.3). The 
top panel shows the reconstructed equation of state w(z ) after 
the PCA analysis. Vertical error bars correspond to mean and 
standard deviations of the q vector parameters, while horizon¬ 
tal error bars are the amplitude of the original binning. The bins 
are not exactly independent but are rather smeared out as illus¬ 
trated in the bottom panel. The bottom panel shows the PCA 
corresponding weights on w{z) as a function of redshift for the 
combination Planck TT+lowP+BSH. In other words, error bars 
in the top panel correspond therefore to the errors in the q param¬ 
eters, which are linear combinations of the p parameters, i.e. a 
smeared out distribution with weights shown in the lower panel. 


matter, as (Huang et al. 2011) 

w = -l + ^e s W— ), (23) 

3 \a de ) 

where the “slope parameter” e s is defined as: 

e s = wl fl =fl de . (24) 

with e v = p/2 being a function of the slope of the po¬ 

tential. Here Mp = 1 / VtVrG is the reduced Planck mass and a de 
is the scale factor where the total matter and DE densities are 
equal. The function F(x) in Eq. (23) is defined as: 


Eq. (23) parameterizes w(a) with one parameter e s , while a de 
depends on Q m and e s and can be derived using an approxi¬ 
mated fitting formula that facilitates numerical computation 
(Huang et al. 2011). Positive (negative) values of e s correspond 
to quintessence (phantom) models. 

Eq. (23) is only valid for late-Universe slow-roll (e v S: 1 
and >] v = M~V "/V <s 1) or the moderate-roll (e v ^ 1 and 
? 7 y < 1) regime. For quintessence models, where the scalar field 
rolls down from a very steep potential, at early times ey(a) » 1, 
however the fractional density £Xi,Ui) —» 0 and the combination 
ey{a)PL^(a) aprroaches a constant, defined to be a second param¬ 
eter = lim^o e v (fl)£Va). 

One could also add a third parameter £ s to capture the time- 
dependence of gy via corrections to the functional dependence 
of w(a) at late time. This parameter is defined as the relative 
difference of d JeyQ^/dy at a — a de and at a —» 0, where y = 

( a/a de ) 312 / Vl + (a/flde) 3 . If too <s 1, 4s is proportional to the 
second derivative of In V(<p), but for large e ra , the dependence is 
more complicated (Huang et al. 2011). In other words, while e s 
is sensitive to the late time evolution of 1 + w(a), e m captures 
its early time behaviour. Quintessence/phantom models can be 
mapped into e s -e m space and the classification can be further 
refined with £ s . For ACDM, all three parameters are zero. 

In Fig. 8 we show the marginalized posterior distribu¬ 
tions at 68.3% and 95.4% confidence levels in the param¬ 
eter space e s -Q m , marginalizing over the other parameters. 
In Fig. 9 we show the current constraints on quintessence 
models projected in e s -6oo space. The constraints are ob¬ 
tained by marginalizing over all other cosmological parameters. 
The models here include exponentials V = Vo exp(-T</>/Mp) 
(Wetterich 1988), cosines from pseudo-Nambu Goldstone 
bosons (pnGB) V = Vo[l + cos(T</>/M P )] (Frieman et al. 
1995; Kaloper & Sorbo 2006), power laws V = Void)lMp) " 
(Ratra & Peebles 1988), and models motivated by supergrav¬ 
ity (SUGRA) V = Vo((/>/Mp)~ a exp [(</>/M P ) 2 ] (Brax & Martin 
1999). The model projection is done with a fiducial D ra = 0.3 
cosmology. We have verified that variations of 1 % compared to 
the fiducial f2 m lead to negligible changes in the constraints. 

Mean values and uncertainties for a selection of cosmo¬ 
logical parameters are shown in Table 2, for both the 1- 
parameter case (i.e., e s only, with e m - 0 and A = 0, de¬ 
scribing “thawing” quintessence/phantom models, where (f> = 
0 in the early Universe) and the 3-parameter case (general 
quintessence/phantom models where an early-Universe fast¬ 
rolling phase is allowed). When we vary the data sets and the¬ 
oretical prior (between the 1-parameter and 3-parameter cases), 
the results are all compatible with ACDM and mutually compat¬ 
ible with each other. Because e s and are correlated, caution 
has to be taken when looking at the marginalized constraints 
in the table. For instance, the constraint on e s is tighter for the 
3-parameter case, because in this case flatter potentials are pre¬ 
ferred in the late Universe in order to slow-down larger cf> from 
the early Universe. A better view of the mutual consistency can 
be obtained from Fig. 9. We find that the addition of polariza¬ 
tion data does not have a large impact on these DE parameters. 
Adding polarization data to Planck+ BSH shifts the mean of e s 
by —1/6 cr and reduces the uncertainty of e s by 20%, while the 
95 % upper bound on e xi remains unchanged. 

5.1.5. Dark energy density at early times 


Fix) = 


Vl + x 3 

X 3/2 


In [xr’F + V1 + x 3 ) 


Quintessence models can be divided into two classes, namely 
(25) cosmologies with or without DE at early times. Although the 
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Parameter Planck+BSH (1-param.) Planck+BSH (3-param.) P/anc/r+WL+BAO/RSD Ptonck+lensing+BSH 


e s . - 0 . 081 °“ - 0 . 111 ° J ° 0 . 141 °'’ - 0 . 03!°; 7 6 

e„. fixed = 0 < 0.76 (95 % CL) < 0.38 (95 % CL) < 0.52 (95 % CL) 

£ s . fixed = 0 not constrained not constrained not constrained 


Table 2. Marginalized mean values and 68 % CL errors for a selection of cosmological parameters for the weakly-coupled scalar 
field parameterization described in the text (Sect. 5.1.4). Here “1-param” in the first column refers to the priors = 0 and £ s = 0 
(slow- or moderate-roll “thawing” models). 


0.6 


0.0 


- 0.6 


- 1.2 


0.24 0.28 0.32 0.36 0.40 

Om 




Fig. 8. Marginalized posterior distributions showing 68 % and 
95 % C.L. constraints on Q m and e s for scalar field models (see 
Sect. 5.1.4). The dashed line for e s = 0 is the ACDM model. 
Here Planck indicates Planck TT+lowP. 


Fig. 10. Marginalized posterior distributions for O e for the early 
DE parameterization in (26) and for different combinations of 
data (see Sect. 5.1.5). Here Planck indicates Planck TT+lowP. 



6qo 


Fig. 9. Marginalized posterior distributions at 68% C.L. and 
95 % C.L. in the parameter space of e s and for scalar field 
modes (see Sect. 5.1.4). We have computed e s and e a , for various 
quintessence potentials V (</;), with the functional forms of V(4>) 
labelled on the figure. The field <p is in reduced Planck mass Mp 
units. The normalization of V{<p) is computed using O m = 0.3. 
Here Planck indicates Planck TT+lowP. 


equation of state and the DE density are related to each other, 
it is often convenient to think directly in terms of DE density 
rather than the equation of state. In this section we provide a 
more direct estimate of how much DE is allowed by the data 
as a function of time. A key parameter for this purpose is O e , 
which measures the amount of DE present at early times (“early 
dark energy,” EDE) (Wetterich 2004). Early DE parameteriza- 
tions encompass features of a large class of dynamical DE mo¬ 
dels. The amount of early DE influences CMB peaks and can be 
strongly constrained when including small-scale measurements 
and CMB lensing. Assuming a constant fraction of f2 e until re¬ 
cent times (Doran & Robbers 2006), the DE density is parame¬ 
terized as: 


OteW = V g o0 3 0) + Qe(1 - fl ~ 3W0) • (26) 

This expression requires two parameters in addition to those of 
ACDM, namely Q e and wo, while O" = 1 - £1'^ is the present 
matter abundance. The strongest constraints to date were dis¬ 
cussed in (Planck Collaboration XVI 2014), finding O e < 0.010 
at 95 % CL using Planck combined with WMAP polarization. 
Here we update the analysis using Planck 2015 data. In Fig. 10 
we show marginalized posterior distributions for Q e for different 
combination of data sets; the corresponding marginalized lim¬ 
its are shown in Table 3, improving substantially current con¬ 
straints, especially when the Planck TT,TE,EE+lowP polariza¬ 
tion is included, leading to Q e < 0.0036 at 95% confidence level 
for Planck TT,TE,EE+lowP+BSH. 
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Parameter Planck TT+lowP Planck TT+lowP Planck TT+lowP Planck TT+lowP Planck TT,TE,EE+lowP 

+lensing+BSH +lensing+WL +lensing+BAO/RSD +lensing+WL+BAO/RSD + BSH 

n c . <0.0071 < 0.0087 <0.0070 <0.0070 <0.0036 

w 0 . <-0.93 <-0.76 <-0.90 <-0.90 <-0.94 


Table 3. Marginalized 95 % limits on O e and wo for the early DE parameterization in Eq. (26) and different combinations of data 
(see Sect. 5.1.5). Including high-7 polarization significantly tightens the bounds. 


As first shown in Pettorino et al. (2013), bounds on O e can be 
weaker if DE is present only over a limited range of redshifts. In 
particular, EDE reduces structure growth in the period after last 
scattering, implying a smaller number of clusters as compared 
to ACDM, and therefore a weaker lensing potential to influence 
the anisotropies at high 7. It is possible to isolate this effect by 
switching on EDE only after last scattering, at a scale factor a e 
(or equivalently for redshifts smaller than z e ). Here we adopt the 
parameterization “EDE3” proposed in Pettorino et al. (2013) to 
which we refer for more details: 


Dde(^) — 


_QleO_ 

QdeO + 3 + G r 0fl 4 

D e 

_QdeO_ 

f^deO 3" O m (|a + Q r 0^ 


for a < a e ; 
for a e < a < a c : 
for a > a c . 


(27) 


In this case, early dark energy is present in the time interval 
a e < a < a c , while outside this interval it behaves as in ACDM, 
including the radiation contribution, unlike in Eq. (26). During 
that interval in time, there is a non-negligible EDE contribu¬ 
tion, parameterized by I2 e . The constant a c is fixed by continuity, 
so that the parameters {f2 e ,a e } fully determine how much EDE 
there was and how long its presence lasted. We choose four fixed 
values of a e corresponding to z e = 10, 50, 200 and 1000 and in¬ 
clude Q e as a free parameter in MCMC runs for each value of 
a e . Results are shown in Fig. 11 where we plot Q e as a function 
of the redshift z e at which DE starts to be non-negligible. The 
smaller the value of z e , the weaker are the constraints, though 
still very tight, with Q e ^ 2 % (95 % CL) for z e ~ 50. 


5.1.6. Compressed likelihood 

Before concluding the set of results on background parame- 
terizations, we discuss here how to reduce the full likelihood 
information to few parameters. As discussed for example in 
Kosowsky et al. (2002) and Wang & Mukherjee (2007), it is pos¬ 
sible to compress a large part of the information contained in 
the CMB power spectrum into just a few numbers 5 : here we 
use specifically the CMB shift parameter R (Efstathiou & Bond 
1999), the angular scale of the sound horizon at last scattering 7 a 
( or equivalently 6, ), as well the baryon density <u h and the scalar 
spectral index n s . The first two quantities are defined as 

R = D a (z*)/c , 7 a = 7rD A (z*)/r s (z t ) = n/6„ , (28) 

where D A (z) is the comoving angular diameter distance to red- 
shift z, z* is the redshift for which the optical depth is unity 
and r s (zt ) = A is the comoving size of the sound horizon at 

5 There are also alternative approaches that compress the power spec¬ 
tra directly, like e.g. PICO (Fendt & Wandelt 2006). 
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Fig. 11. Amount of DE at early times 12,, as a function of the 
redshift z e after which early DE is non-negligible (see Eq. (27), 
Sect. 5.1.5) for different combinations of data sets. The heights 
of the columns give the limit at 95 % CL on Q . e , as obtained from 
Monte Carlo runs for the values z e = 10, 50, 200 and 1000. The 
width of the columns has no physical meaning and is just due to 
plotting purposes. Here Planck indicates Planck TT+lowR 



10.0 100.0 1000.0 


1 + z e 


z*. These numbers are effectively observables and they apply to 
models with either non-zero curvature or a smooth DE compo¬ 
nent (Mukherjee et al. 2008). It should be noted, however, that 
the constraints on these quantities, especially on R, are sensitive 
to changes in the growth of perturbations. This can be seen easi¬ 
ly with the help of the “dark degeneracy” (Kunz 2009), i.e., the 
possibility to absorb part of the dark matter into the dark energy, 
which changes £2 m without affecting observables. For this rea¬ 
son the compressed likelihood presented here cannot be used for 
models with low sound speed or modifications of gravity (and is 
therefore located at the end of this “background” section). 

The marginalized mean values and 68 % confidence inter¬ 
vals for the compressed likelihood values are shown in Table 4 
for Planck TT+lowR The posterior distribution of {R, 7a, a>b, « s } 
is approximately Gaussian, which allows us to specify the like¬ 
lihood easily by giving the mean values and the covariance ma¬ 
trix, as derived from a Monte Carlo Markov chain (MCMC) ap¬ 
proach, in this case from the grid chains for the wCDM model. 
Since these quantities are very close to observables directly 
derivable from the data, and since smoothly parameterized DE 
models are all compatible with the Planck observations to a com¬ 
parable degree, they lead to very similar central values and es¬ 
sentially the same covariance matrix. The Gaussian likelihood 
in {R, 7a, cun, given by Table 4 is thus useful for combin¬ 
ing Planck temperature and low-7 polarization data with other 
data sets and for inclusion in Fisher matrix forecasts for future 
surveys. This is especially useful when interested in parameters 
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Smooth DE models Planck TT+lowP R Cl b h 2 n s 


R . 1.7488 ±0.0074 1.0 0.54 -0.63 -0.86 

C A . 301.76 ±0.14 0.54 1.0 -0.43 -0.48 

Cl b h 2 . 0.02228 ± 0.00023 -0.63 -0.43 1.0 0.58 

n s . 0.9660 ±0.0061 -0.86 -0.48 0.58 1.0 


Marginalized over A L Planck TT+lowP R C A Q b h 2 n s 


R . 1.7382 ±0.0088 1.0 0.64 -0.75 -0.89 

. 301.63 ±0.15 0.64 1.0 -0.55 -0.57 

Cl b h 2 . 0.02262 ± 0.00029 -0.75 -0.55 1.0 0.71 

n s . 0.9741 ±0.0072 -0.89 -0.57 0.71 1.0 


Table 4. Compressed likelihood discussed in Sect. 5.1.6. The left columns give the marginalized mean values and standard deviation 
for the parameters of the compressed likelihood for Planck TT+lowP, while the right columns present the normalized covariance or 
correlation matrix D for the parameters of the compressed likelihood for Planck TT+lowP. The covariance matrix C is then given 
by C ij = CTiCTjDjj (without summation), where <r, is the standard deviation of parameter i. While the upper values were derived for 
n CDVI and are consistent with those of ACDM and the {wo,w a } model, we marginalized over the amplitude of the lensing power 
spectrum for the lower values, which leads to a more conservative compressed likelihood. 


such as {wo,w fl }, for which the posterior is very non-Gaussian 
and cannot be accurately represented by a direct covariance ma¬ 
trix (as can be seen in Fig. 4). 

The quantities that make up the compressed likelihood are 
supposed to be “early universe observables” that describe the ob¬ 
served power spectrum and are insensitive to late time physics. 
However, lensing by large-scale structure has an important 
smoothing effect on the C{ and is detected at over 10 cr in the 
power spectrum (see section 5.2 of Planck Collaboration XIII 
(2015)). We checked by comparing different MCMC chains that 
the compressed likelihood is stable for ACDM, wCDM and the 
{wo, w fl ) model. However, the “geometric degeneracy” in curved 
models is broken significantly by the impact of CMB lensing on 
the power spectrum (see Fig. 25 of Planck Collaboration XIII 
2015) and for non-flat models one needs to be more careful. For 
this reason we also provide the ingredients for the compressed 
likelihood marginalized over the amplitude Al of the lensing 
power spectrum in the lower part of Table 4. Marginalizing 
over Al increases the errors in some variables by over 20 % 
and slightly shifts the mean values, giving a more conserva¬ 
tive choice for models where the impact of CMB lensing on the 
power spectrum is non-negligible. 

We notice that the constraints on {R, a> b , n s ] given in 
Table 4 for Planck TT+lowP data are significantly weaker 
than those predicted by table II of Mukherjee et al. (2008), 
which were based on the “Planck Blue Book” specifications 
Planck Collaboration (2005). This is because these forecasts also 
used high-/ polarization. If we derive the actual Planck covari¬ 
ance matrix for the Planck TT,TE,EE+lowP likelihood then we 
find constraints that are about 50% smaller than those given 
above, and are comparable and even somewhat stronger than 
those quoted in Mukherjee et al. (2008). The mean values have 
of course shifted to represent what Planck has actually mea¬ 
sured. 

5.2. Perturbation parameterizations 

Up to now we have discussed in detail the ensemble of back¬ 
ground parameterizations, in which DE is assumed to be a 
smooth fluid, minimally interacting with gravity. General modi¬ 
fications of gravity, however, change both the background and 
the perturbation equations, allowing for contribution to cluste¬ 
ring (via a sound speed different than unity) and anisotropic 


stress different from zero. Here we illustrate results for pertur¬ 
bation degrees of freedom, approaching MG from two differ¬ 
ent perspectives, as discussed in Sect. 3. First we discuss results 
for EFT cosmologies, with a “top-down” approach that starts 
from the most general action allowed by symmetry and selects 
from there interesting classes belonging to so-called “Horndeski 
models”, which, as mentioned in Sect. 3.2.1, include almost all 
stable scalar-tensor theories, universally coupled, with second- 
order equations of motion in the fields. We then proceed by pa¬ 
rameterizing directly the gravitational potentials and their com¬ 
binations, as illustrated in Sect. 3.2.2. In this way we can test 
more phenomenologically their effect on lensing and clustering, 
in a “bottom-up” approach from observations to theoretical mo¬ 
dels. 


5.2.1. Modified gravity: EFT and Horndeski models 

The first of the two approaches described in Sect. 3.2.1 adopts ef¬ 
fective field theory (EFT) to investigate DE (Gleyzes et al. 2013; 
Gubitosi et al. 2013), based on the action of Eq. (2). The param¬ 
eters that appear in the action, when choosing the nine time- 
dependent functions {O, c,A, ,nv,\, de¬ 

scribe the effective DE. The full background and perturbation 
equations for this action have been implemented in the publicly 
available Boltzmann code EFTCAMB (Hu et al. 2014; Raveri et al. 
2014) 6 . Given an expansion history (which we fix to be ACDM, 
i.e., effectively w= — 1) and an EFT function Q (a), EFTCAMB 
computes c and A from the Friedmann equations and the as¬ 
sumption of spatial flatness (Hu et al. 2014). As we have seen in 
Sect. 5.1, for smooth DE models the constraints on the DE equa¬ 
tion of state are compatible with w — — 1; hence this choice is 
not a limitation for the following analysis. In addition, EFTCAMB 
uses a set of stability criteria in order to specify whether a given 
model is stable and ghost-free, i.e. without negative energy den¬ 
sity for the new degrees of freedom. This will automatically 
place a theoretical prior on the parameter space while perform¬ 
ing the MCMC analysis. 

The remaining six functions, A7 [ \ ;V/j, M\, ;V/j, M 1 , ni 2 r are 
internally redefined in terms of the dimensionless parameters cr, 


6 http://vjww.lorentz.leidenuniv.nl/-hu/codes/, version 
1.1, Oct. 2014. 
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with i running from 1 to 6: 
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We will always demand that 

m\ = 0 (or equivalently cr b = 0), 

M\ = -M? (or equivalently a\ = -af), 


(29) 

(30) 


which eliminates models containing higher-order spatial deriva¬ 
tives (Gleyzes et al. 2013, 2014). In this case the nine functions 
of time discussed above reduce to a minimal set of five functions 
of time that can be labelled {ctm, uk> o' B , aj, an], in addition to 
the Planck mass M 2 (the evolution of which is determined by H 
and a M ), and an additional function of time describing the back¬ 
ground evolution, e.g., H(a). The former are related to the EFT 
functions via the following relations: 


M; = m^Cl + M\\ 

(31) 

M;Ha M = mlCl + M 2 2 ; 

(32) 

M;H 2 a K = 2c + 4M^ 

(33) 

M;Ha B = -/« 2 0 - M,; 

(34) 

M 2 a T = -M 2 ; 

(35) 

M;a H = 2M 2 - m\. 

(36) 

These five a functions are closer to a physical descrip¬ 
tion of the theories under investigation (Bellini & Sawicki 
2014). For example: aj enters in the equation for gravita¬ 
tional waves, affecting their speed and the position of the pri¬ 
mordial peak in /1-mode polarization; aM affects the lensing 
potential, but also the amplitude of the primordial polariza¬ 
tion peak in B-modes (Amendola et al. 2014; Raveri et al. 2014; 
Pettorino & Amendola 2014). It is then possible to relate the de¬ 
sired choice for the Horndeski variables to an appropriate choice 
of the EFT functions, 

d T (Ml) = PlM;a M , 

(37) 

ml(Cl + 1) = (1 + a T )M;, 

(38) 

M\ = -u t M 2 , 

(39) 

4 Ml = M; ( H 1 a^-2c, 

(40) 

M\ = -M;9Ca B + w 2 Q, 

(41) 

2 M~ = M 2 (ar H - aj -), 

(42) 


where 'H is the conformal Hubble function, mo the bare Planck 
mass and M, the effective Plank mass. Fixing % corresponds 
to fixing M* through Eq. (37). Once a t has been chosen, Cl is 
obtained from Eq. (38). Finally, a B determines M 3 via Eq. (41), 
while the choice of an fixes M 2 via Eq. (42). In this way, our 
choice of the EFT functions can be guided by the selection of 
different “physical” scenarios, corresponding to turning on dif¬ 
ferent Horndeski functions. 

To avoid possible consistency issues with higher derivatives, 
we set 7 M 2 — Ml - 0 in order to satisfy Eq. (30). From Eq. (39) 
and Eq. (31) this implies aj = 0, so that tensor waves move with 

7 Because of the way EFTCAMB currently implements these equations 
internally, it is not possible to satisfy Eq. (30) otherwise. 


the speed of light. In addition, we set an = 0 so as to remain 
in the original class of Horndeski theories. As a consequence, 
M 2 = 0 from Eq. (42) and M 2 = m^( 1 + Cl) from Eq. (31). For 
simplicity we also turn off all other higher-order EFT operators 
and set M 3 = M\ = 0. Comparing Eq. (32) and Eq. (41), this 
implies a B = -Q"m- 

In summary, in the following we consider Horndeski models 
in which q-m = -ctb, a k is fixed by Eq. (34), with Mi = 0 as a 
function of c and aj = an = 0. We are thus considering non- 
minimally coupled “K-essence” type models, similar to the ones 
discussed in Sawicki et al. (2013). 

The only free function in this case is o'm, which is linked to 
Q through: 


O'M 


a d£2 
£2 + 1 da 


(43) 


By choosing a non-zero um (and therefore a time evolving Cl) 
we introduce a non-minimal coupling in the action (see Eq. 2), 
which will lead to non-zero anisotropic stress and to modifica¬ 
tions of the lensing potential, typical signatures of MG models. 
Here we will use a scaling ansatz, o'm = cimoo! 3 , where ctmo is 
the value of <jm today, and [1 > 0 determines how quickly the 
modification of gravity decreases in the past. 

Integrating Eq. (43) we obtain 


Cl(a) = exp | ^ — 1, (44) 

which coincides with the built-in exponential model of EFTCAMB 
for I>o = o' mo / [i. The marginalized posterior distributions for 
the two parameters Clo and ji are plotted in Fig. 12 for different 
combinations of data. For umo = 0 we recover ACDM. For small 
values of £2 q and for/? = 1, the exponential reduces to the built-in 
linear evolution in EFTCAMB, 


Cl(a) = Clo a. (45) 

The results of the MCMC analysis are shown in Table 5. For 
both the exponential and the linear model we use a flat prior 
Cl 0 6 [0-1] - For the scaling exponent /? of the exponential model 
we use a flat prior /? e (0,3]. For /? —» 0 the MG parameter aM 
remains constant and does not go to zero in the early Universe, 
while for /? = 3 the scaling would correspond to M functions in 
the action (2) which are of the same order as the relative energy 
density between DE and the dark matter background, similar to 
the suggestion in Bellini & Sawicki (2014). An important fea¬ 
ture visible in Fig. 12 is the sharp cutoff at /? « 1.5. This cutoff 
is due to “viability conditions” that are enforced by EFTCAMB 
and that reject models due to a set of theoretical criteria (see 
Hu et al. (2014) for a full list of theoretical priors implemented 
in EFTCAMB). Disabling some of these conditions allows to ex¬ 
tend the acceptable model space to larger /?, and we find that the 
constraints on a mo continue to weaken as /? grows further, ex¬ 
tending Fig. 12 in the obvious way. We prefer however to use 
here the current public EFTCAMB version without modifications. 
A better understanding of whether all stability conditions imple¬ 
mented in the code are really necessary or exclude a larger region 
than necessary in parameter space will have to be addressed in 
the future. The posterior distribution of the linear evolution for Cl 
is shown in Fig. 13 and is compatible with ACDM. Finally, it is 
interesting to note that in both the exponential and the linear ex¬ 
pansion, the inclusion of WL data set weakens constraints with 
respect to Planck TT+lowP alone. This is due to the fact that 
in these EFT theories, WL and Planck TT+lowP are in tension 
with each other, WL preferring higher values of the expansion 
rate with respect to Planck. 
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Parameter TT+lowP+BSH TT+lowP+WL TT+lowP+BAO/RSD TT+lowP+BAO/RSD+WL TT,TE,EE TT,TE,EE+BSH 


Linear EFT. 

ff M o. < 0.052(95 %CL) < 0.072(95 %CL) < 0.057(95 %CL) < 0.074(95 %CL) < 0.050(95 %CL) < 0.043(95 %CL) 

H 0 . 67.69 ±0.55 67.75 ±0.95 67.63 ± 0.63 67.89 ±0.62 67.17 ±0.66 67.60 ±0.48 

erg. 0.826 ±0.015 0.818 ±0.014 0.822 ±0.014 0.814 ±0.014 0.830 ±0.013 0.830 ±0.014 


Exponential EFT . 

q-mo . < 0.063(95 %CL) < 0.092(95 %CL) < 0.066(95 %CL) < 0.097(95 %CL) < 0.054(95 %CL) < 0.062(95 %CL) 

/ ; . o- 87 ^ o,)2 

H 0 . 67.70 ±0.56 67.78 ±0.96 67.60 ±0.62 67.87 ± 0.63 67.15 ±0.65 67.58 ± 0.46 

cr 8 . 0.826 ±0.015 0.817 ±0.014 0.821 ± 0.014 0.814 ±0.014 0.830 ±0.013 0.830 ±0.013 


Table 5. Marginalized mean values and 68 % CL intervals for the EFT parameters, both in the linear model, ctmo, and in the 
exponential one, {auo,P} (see Sect. 5.2.1). Adding CMB lensing does not improve the constraints, while small-scale polarization 
can more strongly constraint aMO- 



«M0 


Fig. 12. Marginalized posterior distributions at 68 % and 95 % 
C.L. for the two parameters a M0 and [5 of the exponential evo¬ 
lution, f 2(a) = exp(Qo^) - 1.0, see Sect. 5.2.1. Here otmo is 
defined as Qo/j and the background is fixed to ACDM. 0\i« = 0 
corresponds to the ACDM model also at perturbation level. Note 
that Planck means Planck TT+lowP. Adding WL to the data sets 
results in broader contours, as a consequence of the slight ten¬ 
sion between the Planck and WL data sets. 



Fig. 13. Marginalized posterior distribution of the linear EFT 
model background parameter, Q, with Q parameterized as a lin¬ 
ear function of the scale factor, i.e., O (a) = or mo a, see Sect. 
5.2.1. The equation of state parameter Wd e is fixed to -1, and 
therefore, Do = 0 will correspond to the ACDM model. Here 
Planck means Planck TT+lowP. Adding CMB lensing to the 
data sets does not change the results significantly; high-f po¬ 
larization tightens the constraints by a few percent, as shown in 
Tab. 5. 


5.2.2. Modified gravity and the gravitational potentials 

The second approach used in this paper to address MG is more 
phenomenological and, as described in Sect. 3.2.2, starts from 
directly parameterizing the functions of the gravitational poten¬ 
tials listed in Eqs. (3)—(6). Any choice of two of those functions 
will fully parameterize the deviations of the perturbations from 
a smooth DE model and describe the cosmological observables 
of an MG model. 

In Simpson et al. (2013) the amplitude of the deviation with 
respect to ACDM was parameterized similarly to the DE-related 
case that we will define as case 1 below, but using /a(a) and 
2(a) instead of p(a,k) and //(a, k f. They found the constraints 
p 0 - 1 = 0.05 ± 0.25 and S 0 - 1 = 0.00 ± 0.14 using RSD 
data from the WiggleZ Dark Energy Survey (Blake et al. 201 1) 


x The parameterization of p and E in Simpson et al. (2013) uses 
Dde(°)/Dde instead of £2 DE (a); theirp 0 and £o correspond to our p 0 - 1 
and Eo - 1 respectively. 


and 6dF Galaxy Survey (6dFGS) (Beutler et al. 2012), together 
with CFHTLenS WL. Baker et al. (2014b) provided forecasts on 
Ho - 1 and E () - 1 for a future experiment that combines galaxy 
clustering and tomographic weak lensing measurements. The 
amplitude of departures from the standard values was parameter¬ 
ized as in Simpson et al. (2013), but a possible scale dependence 
was introduced. In Zhao et al. (2010), the authors constrained ho 
and 7/0 and derived from those the limits on Eo, using WMAP-5 
data along with CFHTLenS and ISW data. Together with a prin¬ 
cipal component analysis, they also constrained p and // assum¬ 
ing a time evolution of the two functions, introducing a transi¬ 
tion redshift z s where the functions move smoothly from an early 
time value to a late time one; they obtained ho = 1.1^'®, t/o = 
0.98+°; 73 forz s = 1 and p 0 = 0.87 ± 0.12, t/ 0 = 1.3 + 0.35 for 
z s = 2. A similar parametrization was also used in Daniel et al. 
( 2010 ) in terms of ho and m (equivalent to ho ~ 1 and 770 -1 in our 
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convention) using WMAP5, Union2, COSMOS and CFHTLenS 
data, both binning these functions in redshift and assuming a 
time evolution (different from the one we will assume in the fol¬ 
lowing), obtaining -0.83 < no <2.1 and -1.6 < m < 2.7 at 
95% confidence level for their present values. In Macaulay et al. 
(2013) the authors instead parameterized 'P/<P (the inverse of 
77 ) as (1 - £) and use RSD data from 6 dFGS, BOSS, LRG, 
WiggleZ and VIPERS galaxy redshift surveys to constrain de¬ 
partures from ACDM; they did not assume a functional form for 
the time evolution of but rather constrained its value at two 
different redshifts (z = 0 and z = 1 ), finding a 2 cr tension with 
the ACDM limit (£ = 0) at z = 1. 

In this paper, we choose the pair of functions n(a, k) (related 
to the Poisson equation for V F) and i](a, k) (related to the gravita¬ 
tional slip), as defined in Eqs. (4) and ( 6 ), since these are the 
functions directly implemented in the publicly available code 
MGCAMB 9 (Zhao et al. 2009; Hojjati et al. 2011) integrated in the 
latest version of CosmoMC. 

Other functional choices can be easily derived from them 
(Baker et al. 2014b). We then parameterize 7 / and // as follows. 
Since the Planck CMB data span three orders of magnitude in (, 
it seems sensible to allow for two scales to be present: 


H(a,k) = 1 + fi(a) 


rj(a, k) = 1 + f 2 (a ) 


1 + c\(AH/k ) 2 
1 + ( AH/k) 2 ’ 
1 + c 2 (AH/k ) 2 
1 + (AH/k) 2 


(46) 

(47) 


For large length scales (small k ), the two functions reduce to 
7/ —* 1 + f\{a)c\ and 77 —> 1 + f 2 (a)c 2 \ for small length scales 
(large k ), one has /j —> 1 +f\(a) and 77 —> 1 +f 2 {a). In other words, 
we implement scale dependence in a minimal way, allowing 77 
and 77 to go to two different limits for small and large scales. Here 
the f are functions of time only, while the q and A parameters 
are constants. The c, give us information on the scale dependence 
of 7 / and 77 , but the f measure the amplitude of the deviation from 
standard GR, corresponding to /j = 77 = 1. 

We choose to parameterize the time dependence of the /j(a) 
functions as 


1. coefficients related to the DE density, f(a) = 

2 . time-related evolution, f(a ) = An + £ 12(1 - a). 

The first choice is motivated by the expectation that the contri¬ 
bution of MG to clustering and to the anisotropic stress is pro¬ 
portional to its effective energy density, as is the case for matter 
and relativistic particles. The second parameterization provides 
a complementary approach to the first: £;i describes the MG 
contribution at late times, while £; 2 is relevant at early times. 
Therefore the adoption of the time-related evolution allows, in 
principle, for deviations from the standard behaviour also at high 
redshift, while the parameterization connected to the DE density 
leads by definition to ( 77 , 77 ) —> 1 at high redshift, since the red¬ 
shift evolution is tied to that of Gde(-c)- 

For case 1 (referred to as “DE-related” parameterization) we 
then have five free parameters, £ 11 , c\, E 22 , c 2 , and A, while 
for case 2 (the “time-related” parameterization) we have two 
additional parameters, £12 and £ 21 . The choice above looks 
very similar to the BZ parameterization (Bertschinger & Zukin 
2008) for the quasi-static limit of f(R) and scalar-tensor theo¬ 
ries. However, we emphasise that Eqs. (46)—(47) should not be 

9 Available at http://www.sfu.ca/~aha25/MGCAHB.html (Feb. 
2014 version), see appendix A of (Zhao et al. 2009) for a detailed de¬ 
scription of the implementation. 


seen as a quasi-static limit of any specific theory, but rather as 
a (minimal) way to allow for (arbitrary) scale dependence, since 
the data cover a sufficiently wide range of scales. Analogously 
to the EFT approach discussed in the previous section, we set 
the background evolution to be the same as in ACDM, so that 
w — -1. In this way the additional parameters purely probe the 
perturbations. 

The effect of the £„ parameters on the CMB temperature and 
lensing potential power spectra has been shown in Fig. 1 for the 
“DE-related” choice. In the temperature spectrum the amplitude 
of the ISW effect is modified; the lensing potential changes more 
than the temperature spectrum for the same amplitude of the E n 
parameter. 

We ran Monte Carlo simulations to compare the theoretical 
predictions with different combinations of the data for both cases 
1 and 2. For both choices we tested whether scale dependence 
plays a role (via the parameters c and A) with respect to the 
scale-independent case in which we fix ci - c 2 = 1. Results 
show that a scale dependence of 7 / and 77 does not lead to a 
significantly smaller x 2 with respect to the scale-independent 
case, both for the DE-related and time-related parameterizations. 
Therefore there is no gain in adding c, and A as extra degrees of 
freedom. For this reason, in the following we will mainly show 
results obtained for the scale-independent parameterization. 

Table 6 shows results for the DE-related case for different 
combinations of the data. Adding the BSH data sets to the Planck 
TT+lowP data does not significantly increase the constraining 
power on MG parameters; Planck polarization also has little im¬ 
pact. On the contrary, the addition of RSD data tightens the con¬ 
straints significantly. The WL contours, including the ultracon¬ 
servative cut that removes dependence on nonlinear physics, re¬ 
sult in weaker constraints. In the table, no — 1 and 770 - 1 are 
obtained by reconstructing Eqs. (46) and (47) from E\\ and E 22 
at the present time. In addition, the present value of the £ pa¬ 
rameter, defined in Eq. (5), can be obtained from 77 and 77 as 
£ = (///2)(1 + 77 ) using Eqs. (4) and ( 6 ). 


DE-related 

- 1 - 1 - 1 - 

Planck 

1.6 - 

1 Planck+VJL 



0.00 0.25 0.50 0.75 1.00 

Eo-1 


Fig. 15. Marginalized posterior distributions for 68 % and 95% 
C.L. for the two parameters {/io — 1, £0 — 1} obtained by evalu¬ 
ating Eqs. (46) and (47) at the present time in the DE-related 
parametrization when no scale dependence is considered (see 
Sect. 5.2.2). £ is obtained as £ = (/r/2)(l + 77 ). The time-related 
evolution would give similar contours. In the labels, Planck 
stands for Planck TT+lowR 
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Parameter 

Planck TT±lowP Planck TT±lowP 

Planck TT±lowP 

Planck TT±lowP 

Planck TT±lowP 

Planck TT,TE,EE+lowP 


±BSH 

±WL 

± BAO/RSD 

±WL±BAO/RSD 

±BSH 


E 11. 0.099$“ 0.06$ 32 

E 22 . 0.99 ±1.3 1.03 ±1.3 

P»- I. 0.07$“ «.04'!S 

770 - 1 . 0.70 ± 0.94 0.72 ± 0.90 

X„ - I . 0.28 ±0.15 0.27 ±0.14 

r . 0.065 ±0.021 0.063 ± 0.020 

Mq (km/s/Mpc) . 68.5 ±1.1 68.17 ±0.58 

c-8 . 0.817$““ 0.816$““ 


-0 20 + ° 19 
u - iu -0.47 

-0 24 +019 
u - -0.33 

1 

O 

U) 

0 

1 + 

© p 

O OO 

1 92 +L4 

^ -0.96 

1.77 ±0.88 

2.07 ± 0.85 

-0 14+° 13 

U - 1 -0.34 

-0 17+ 014 
'-0.23 

-0 21 +012 
0-4.1-021 

1 36 +lo 

• -0.69 

1.23 ±0.62 

1.45 ±0.60 

0 34 +017 

-0.14 

0.29 ±0.13 

0.31 ±0.13 

0 061 + ° 020 
u.uui 00 22 

0.062 ±0.019 

0.057 ± 0.019 

69.2 ± 1.1 

68.26 ± 0.69 

68.55 ± 0.66 

0-786^ 2 ‘ 

O 7Q9+0.021 

u. 1 y-_ og25 

0-781/^023 


0 - 08 !°;” 

0.9 ± 1.2 

0 06 +0 ' 23 
u - uu -0.48 

0.60 ±0.86 

0.23 ±0.13 

0.060 ±0.019 

67.90 ± 0.48 

0 - 816 t”[ 


Table 6. Marginalized mean values and 68 % C.L. errors on cosmological parameters and the parameterizations of Eqs. (46) and 
(47) in the DE-related case (see Sect. 5.2.2), for the scale-independent case. 


Max. degeneracy 

Planck TT±lowP 

Planck TT±lowP 

Planck TT±lowP 

Planck TT±lowP 

Planck TT±lowP 




±BSH 

±WL 

±BAO/RSD 

±WL±BAO/RSD 

DE-related .... 

0.84/«; 3 ° 

(2.1(7) 

0-80/q39 (2.1*7-) 

1-08 4 q42 (2.6*7) 

0.90$ 33 (2.4 0 -) 

1.03 ±0.34 (3.0(7) 

+ CMB lensing 

0.42Z 9 ;* 3 

(1.2*7-) 

0.38$ 3 * (1.4*7-) 

0.58$“ (1.6*7) 

0.40$“ (1.4*7) 

0 51 +0 - 21 

U.J1-Q30 

(1.7(7) 

Time-related . . . 

0-67$“ 

(1.0*7-) 

0-69/^67 (1.0*7-) 

1.12$S ( 1 - 8o ‘) 

0.55$ 23 (1.7(7) 

0 70 +0 - 27 
yj ‘' - 0.33 

(2.1(7) 


Table 7. Marginalized mean values and 68% C.L. errors on the present day value of the function 2 \ju(z,k) - 1] + [r/(z,k) - 1], 
which corresponds to the (approximate) maximum degeneracy line identified within the 2 dimensional posterior distributions. This 
function gives a quick idea of the maximum possible tension found for each data set combination in these classes of models, for the 
scale-independent case. The upper part of the table refers to the DE-related parametrisation, with and without CMB lensing, while 
the lower part refers to the time-related one (see Sect. 5.2.2). For convenience, we write explicitly in brackets for each case the 
tension in units of cr with respect to the standard ACDM zero value. The DE-related case is more in tension than the time-related 
parameterization, with a maximum tension that ranges between 2.1 cr and 3cr, depending on the data sets. When CMB lensing is 
included, also the DE-related parameterization becomes compatible with ACDM, with a maximum possible ‘tension’ of at most 1.7 
cr when WL and BAO/RSD are included. 
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Fig. 14. 68 % and 95 % contour plots for the two parameters j/io - 1 ,770 - 1) obtained by evaluating Eqs. (46) and (47) at the present 
time when no scale dependence is considered (see Sect. 5.2.2). We consider both the DE-related (left panel) and time-related 
evolution cases (right panel). Results are shown for the scale-independent case (ci = - 1). In the labels, Planck stands for Planck 

TT±lowP. 


Some tension appears, in particular, when plotting the 
marginalized posterior distributions in the planes (yUo - 1 ,770 - 1) 
and (710 - 1, X ( ) - 1), as shown in Figs. 14 and 15. Here the con¬ 


straints on the two parameters that describe the perturbations in 
MG are simultaneously taken into account. In Fig. 14, left and 
right panels refer to the DE-related and time-related parameteri- 
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DE-related 



Time-related 


Planck+BSH 

P/anc/c+WL+BAO/RSD 



Fig. 16. Redshift dependence of the function 2 [p(z, k) - 1] + [r](z, k) - 1], defined in Eqs. (46,47), which corresponds to the maximum 
degeneracy line identified within the 2 dimensional posterior distributions. This combination shows the strongest allowed tension 
with ACDM. The left panel refers to the DE-related case while the right panel refers to the time-related evolution (see Sect. 5.2.2). 
In both panels, no scale dependence is considered. The coloured areas show the regions containing 68 % and 95 % of the models. In 
the labels, Planck stands for Planck TT+lowP. 
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Fig. 17. 68 % and 95 % marginalised posterior distributions for the two parameters {/Jo - 1, 770 - 1) obtained by evaluating Eqs. (46) 
and (47) at the present time when no scale dependence is considered (see Sect. 5.2.2). Here we show the effect of CMB lensing, 
which shifts the contours towards ACDM. In the labels, Planck stands for Planck TT+lowP. 


zations defined in 5.2.2, respectively, while the dashed lines in¬ 
dicate the values predicted in ACDM. Interestingly, results ap¬ 
pear similar in both parameterizations. In the DE-related case 
(left panel), the ACDM point lies at the border of the 2 cr con¬ 
tour, already when considering Planck TT+lowP alone. More 
precisely, when looking at the goodness of fit, with respect to the 
standard ACDM assumption, the MG scenario (which includes 
two extra parameters An and £ 22 ) leads to an improvement of 
A+ = -6.3 when using Planck TT+lowP (similarly divided be¬ 
tween lowP and TT) and of Ax 2 = -6.4 when including BSH 

(with a Axq MB -5.6 equally divided between TT and lowP). 

When Planck data (TT+lowP) are combined also with WL data, 
the tension increases to Ax 2 = -10.6 (with the CMB still con¬ 
tributing about the same amount, A+ MB = -6.0). When con¬ 
sidering Planck TT+lowP+BAO/RSD, Ay 2 = -8.1 with respect 


to ACDM while, when combining both WL and BAO/RSD, the 
tension is maximal, with A+ = -10.8 and Xqmb = ”6.9. There 
is instead less tension for the time-related parameterization, as is 
visible in the right panel of Fig. 14. 

Once the behaviour of the coefficients in the two param¬ 
eterizations is known, we can use Eq. (46) to reconstruct the 
evolution of /; and // with scale factor (or redshift, equiva¬ 
lently). In Fig. 16 we choose to show the linear combination 
2[fj(z, k) - 1] + [rj(z, k) - 1], which corresponds approximately to 
the maximum degeneracy line in the 2 dimensional /j - 1, 77 - 1 
parameter space, which allows to better visualize the joint con¬ 
straints on /; and // and their maximal allowed departure from 
ACDM. As expected, the DE-related dependence forces the 
combination to be compatible with ACDM in the past, when the 
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DE density is negligible; the time-related parameterization, in¬ 
stead allows for a larger variation in the past. 

The tension can be understood by noticing that the best fit 
power spectrum corresponds to a value of p and rj (P\\ = -0.3, 
E 22 -2.2 for Planck TT+lowP) close to the thick long dashed 
line shown in Fig. 1 for demonstration. This model leads to 
less power in the CMB at large scales and a higher lensing 
potential, which is slightly preferred by the data points with 
respect to ACDM. This explains also why the MG parame¬ 
ters are somewhat degenerate with the lensing amplitude Al 
( which is an ‘unphysical’ parameter redefining the lensing am¬ 
plitude that affects the CMB power spectrum). As discussed 
in Planck Collaboration XIII (2015) (see for ex. Sect. 5.1.2), 
ACDM would lead to a value of Al (Calabrese et al. 2008) 
somewhat larger than 1. When varying it in MG, we find a 
mean value of Al = 1.116^gj? 5 which is compatible with 
Al = 1 at 1 cr. The price to pay is the tension with ACDM in 
MG parameter space, which compensates the need for a higher 
Al that one would have in ACDM. The CMB lensing likeli¬ 
hood extracted from the 4- point function of the Planck maps 
(Planck Collaboration XV 2015) on the other hand does not pre¬ 
fer a higher lensing potential and agrees well with ACDM. For 
this reason the tension is reduced when we add CMB lensing, 
as shown in Fig. 17. We also note that constraints for this class 
of model are sensitive to the estimation of the optical depth r. 
Smaller values of t tend to shift the results further away from 
ACDM. 

In order to have a quick overall estimate of the tension 
for all cases discussed above, we then show in Table 7 the 
marginalized mean and 68% CF errors for the linear combina¬ 
tion 2[p(z, k) - 1] + [rj(z, k) - 1] that . In the table, we indicate 
in brackets, for convenience, the ‘tension’ with ACDM for each 
case. This is the maximum allowed tension, since it is calcu¬ 
lated along the maximum degeneracy direction. The DE-related 
parameterization is more in tension with ACDM than the time- 
related one. The maximum tension reaches 3 cr when including 
WF and BAO/RSD, being therefore mainly driven by external 
data sets. The inclusion of CMB lensing shifts results towards 
ACDM, as discussed. 

Finally, in general, p and // depend not only on redshift but 
also on scale, via the parameters (c,. A). When marginalizing 
over them, constraints become weaker, as expected. The com¬ 
parison with the scale-independent case is shown in Fig. 18 for 
Planck TT+lowP+BSH and different values of k. When allow¬ 
ing for scale dependence, the tension with ACDM is washed out 
by the weakening of the constraints and the goodness of fit does 
not improve with respect to the scale independent case. 

5.3. Further examples of particular models 

Quite generally, DE and MG theories deal with at least one extra 
degree of freedom that can usually be associated with a scalar 
field. For ‘standard’ DE theories the scalar field couples mini¬ 
mally to gravity, while in MG theories the field can be seen as the 
mediator of a fifth force in addition to standard interactions. This 
happens in scalar-tensor theories (including f(R) cosmologies), 
massive gravity, and all coupled DE models, both when matter 
is involved or when neutrino evolution is affected. Interactions 
and fifth forces are therefore a common characteristic of many 
proposed models, the difference being whether the interaction is 
universal (i.e., affecting all species with the same coupling, as in 
scalar-tensor theories) or is different for each species (as in cou¬ 
pled DE, Wetterich 1995b; Amendola 2000 or growing neutrino 
models, Fardon et al. 2004; Amendola et al. 2008a). In the fol¬ 


lowing we will test well known examples of particular models 
within all these classes. 

5.3.1. Minimally coupled DE: sound speed and k-essence 

In minimally coupled quintessence models, the sound speed 
is cl — 1 and DE does not contribute significantly to clus¬ 
tering. However, in so-called “k-essence” models, the kinetic 
term in the action is generalised to an arbitrary function of 
(V0) 2 (Armendariz-Picon et al. 2000): the sound speed can then 
be different from the speed of light and if c s <s 1, the DE per¬ 
turbations can become non-negligible on sub-horizon scales and 
impact structure formation. To test this scenario we have per¬ 
formed a series of analyses where we allow for a constant equa¬ 
tion of state parameter w and a constant speed of sound c 2 (with 
a uniform prior in logc s ). We find that the limits on w do not 
change from the quintessence case and that there is no signifi¬ 
cant constraint on the DE speed of sound using current data. This 
can be understood as follows: on scales larger than the sound 
horizon and for w close to -1, DE perturbations are related to 
dark matter perturbations through Ade — (1 + w) A m /4 and in¬ 
side the sound horizon they stop growing because of pressure 
support (see e.g., Creminelli et al. 2009; Sapone & Kunz 2009). 
In addition, at early times the DE density is much smaller than 
the matter density, with pDE/Pm = [(1 - D m )/n m ]Gr 3w '. Since 
the relative DE contribution to the perturbation variable Q(a, k) 
defined in Eq. (3) scales like pDEADE/(PmA m ), in k-essence type 
models the impact of the DE perturbations on the total clustering 
is small when 1 + w « 0. For the DE perturbations in k-essence to 
be detectable, the sound speed would have had to be very small, 
and |1 + w\ relatively large. 

5.3.2. Massive gravity and generalized scalar field models 

We now give two examples of subclasses of Horndeski models, 
written in terms of an alternative pair of DE perturbation func- 
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Fig. 18. 68 % and 95 % contour plots for the two parameters 
{po(k) - l,rjo(k) - 1) obtained by evaluating Eqs. (46) and 
(47) at the present time for the DE-related parameterization 
(see Sect. 5.2.2). We consider both the scale-independent and 
scale-dependent cases, choosing k values of 10 _10 Mpc 1 and 
10 2 Mpc _1 . 
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tions (with respect to p and 77 used before, for example), given 
by the anisotropic stress cr and the entropy perturbation T: 

(48 ) 

P dp 

When F = 0 the perturbations are adiabatic, that is 6p = 

For this purpose, it is convenient to adopt the ‘equation of 
state’ approach described in Battye et al. (2015); Soergel et al. 
(2015). The gauge-invariant quantities T and cr can be specified 
in terms of the other perturbation variables, namely 6p, 9, h and 
77 in the scalar sector, and their derivatives. 

We then show results for two limiting cases in this for¬ 
malism, corresponding to Lorentz-violating massive gravity 
(LVMG) for which (cr ± 0,T = 0) and generalized scalar field 
models (GSF) in which the anisotropic stress is zero (cr = 0, T V 
0). 


This has three extra parameters (a,/3\,/32), in addition to w. If 
f5\ - \ and [h = 0 this becomes the generalized k-essence 
model. An example of this class of models is “kinetic gravity 
braiding” (Deffayet et al. 2010) and similar to the non-minimally 
coupled k-essence discussed via EFT in Sect. 5.2.1. The a pa¬ 
rameter in Eq. (50) can be now interpreted as a sound speed, 
unconstrained as in results above. There are however two addi¬ 
tional degrees of freedom, [3\ and /T. RSD data are able to con¬ 
strain them, with the addition of Planck lensing and WL making 
only a minor change to the joint constraints. As in the LVMG 
case, we use a new basis y,- = 11 + w\ a ‘[3j in the MCMC analysis, 
where a\ = 0.2 and o-o = 1 were chosen to decorrelate w and 
y,-. The resulting 2 cr upper limits are yi < 0.67 and yi < 0.61 
(for w > -1), J 2 < 2.4 (for w < -1) for the combination of 
Planck TT+lowP+lensing+BAO/RSD+WL. As for the LVMG 
case, there is no detection of a deviation from ACDM and for 
w — -1 there are no constraints on [j\ and /T. 


Lorentz-violating massive gravity (LVMG) If the Lagrangian is 
£. = -Cigjif ) (i.e. only written in terms of metric perturbations, 
as in the EFT action) and one imposes time translation invari¬ 
ance (but not spatial translational invariance), one finds that this 
corresponds to an extra degree of freedom, , that has a physi¬ 
cal interpretation as an elastic medium, or as Lorentz-violating 
massive gravity (Dubovsky 2004; Rubakov & Tinyakov 2008; 
Battye & Pearson 2013). In this case, the scalar equations are 
characterized by T = 0 (the model is adiabatic) and a non¬ 
vanishing anisotropic stress: 


cr = (w- c 2 ) 


6 

1 + w 



(49) 


including one degree of freedom, the sound speed c 2 , which 
can be related equivalently to the shear modulus of the elas¬ 
tic medium or the Lorentz violating mass. Tensor (gravitational 
wave) equations will also include a mass term. The low sound 
speed may lead to clustering of the DE fluid, which allows the 
data to place constraints on c 2 . But as w approaches -1, the DE 
perturbations are suppressed and the limits on the sound speed 
weaken. We can take this degeneracy between 1 +w and c 2 into 
account by using the combination A c = |1 + wriog^c;: in the 
MCMC analysis, where a = 0.35 was chosen to decorrelate 
w and A c . With this, we find Planck TT+lowP+lensing gives 
lower limit of A c > -1.6 at 2<x and a tighter one when includ¬ 
ing BAO/RSD and WL, with A c > -1.3 at 2cr. For any w 4 -1 
these limits can be translated into limits on log 10 c 2 by comput¬ 
ing Ad\\ + w\ a . The ACDM limit is however fully compatible 
with the data, i.e. there is no detection of any deviation from 
w - — 1 (and in this limit c 2 is unconstrained). 


Generalized scalar field models (GSF) One can allow for 
generalized scalar fields by considering a Lagrangian X = 

in which the dependence on the 
scalar fields is made explicit, imposing full reparameteriza¬ 
tion invariance (x 1 ' —> x v + X 1 ), allowing for only linear cou¬ 
plings in d a g MV and second-order field equations. In this case the 
anisotropic stresses are zero and 


wF = (a — w) 


6-3/3iX((l + w)9- 


3PXH(\ + w) ■ 
2k 2 - 6 (<H - Xt 2 ) 


3(1-yS 2 -/jQTm+w) .A 

6 XI - 18'H'W + &H 2 + 2 k 2 X( 7 ' 


5.3.3. Universal couplings: f(R) cosmologies 


A well-investigated class of MG models is constituted by the 
f(R) theories that modify the Einstein-Hilbert action by substi¬ 
tuting the Ricci scalar with a more general function of itself: 

5 " + /(*)) + f d 4 xX M (x l ,g l nX (51) 

where k 2 = 87 rG. f(R) cosmologies can be mapped to a subclass 
of scalar-tensor theories, where the coupling of the scalar field 
to the matter fields is universal. 

For a fixed background, the Friedmann equation provides 
a second-order differential equation for f(R(a)) (see e.g., 
Song et al. 2007a; Pogosian & Silvestri 2008). One of the initial 
conditions is usually set by requiring 


lim 

/?—» OO 


m 

R 


= 0, 


(52) 


and the other initial (or boundary condition), usually called B 0 , 
is the present value of 


B{z) = 


fm HR 
1 + f R H - H 2 ' 


(53) 


Here, /« and /«* are the first and second derivatives of f(R), 
and a dot means a derivative with respect to conformal time. 
Higher values of Bq suppress power at large scales in the 
CMB power spectrum, due to a change in the ISW effect. This 
also changes the CMB lensing potential, resulting in slightly 
smoother peaks at higher £s (Song et al. 2007b; Schmidt 2008; 
Bertschinger & Zukin 2008; Marchini et al. 2013). 

It is possible to restrict EFTcamb to describe f(R)- 
cosmologies. Given an evolution history for the scale factor and 
the value of Bo, EFTcamb effectively solves the Friedmann equa¬ 
tion for f(R). It then uses this function at the perturbation level 
to evolve the metric potentials and matter fields. The merit of 
EFTcamb over the other available similar codes is that it checks 
the model against some stability criteria and does not assume the 
quasi-static regime, where the scales of interest are still linear but 
smaller than the horizon and the time derivatives are ignored. 

As shown in Fig. 19, there is a degeneracy between the op¬ 
tical depth, r, and the f(R) parameter. Bo- Adding any structure 
formation probe, such as WL, RSD or CMB lensing, breaks the 
degeneracy. Figure 20 shows the likelihood of the B t) param¬ 
eter using EFTcamb, where a ACDM background evolution is 
assumed, i.e., wde = — 1 . 
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As the different data sets provide constraints on Bo that 
vary by more than four orders of magnitude, we show plots for 
log 10 B()\ to make these figures we use a uniform prior in log 10 Bq 
to avoid distorting the posterior due to prior effects. However, 
for the limits quoted in the tables we use Bq (without log) as the 
fundamental quantity and quote 95 % limits based on Bo. In this 
way the upper limit on Bo is effectively given by the location of 
the drop in probability visible in the figures, but not influenced 
by the choice of a lower limit of log 10 By. Overall this appears 
to be the best compromise to present the constraints on the Bq 
parameter. In the plots, the GR value (By = 0) is reached by a 
plateau stretching towards minus infinity. 



-7.5 -6.0 -4.5 -3.0 -1.5 0.0 


log(B 0 ) 

Fig. 19. 68 % and 95 % contour plots for the two parameters, 
|Log I 0 (By), r) (see Sect. 5.3.3). There is a degeneracy between 
the two parameters for Planck TT+lowP+BSH. Adding lensing 
will break the degeneracy between the two. Here Planck indi¬ 
cates Planck TT+lowR 


Finally, we note that f(R ) models can be studied also with 
the MGcamb parametrization, assuming the quasi-static limit. We 
find that for the allowed range of the Bo parameter, the re¬ 
sults with and without the quasi-static approximation are the 
same within the uncertainties. The 95 % confidence intervals 
are reported in Table 8 . These values show an improvement 
over the WMAP analysis made with MGcamb (Bo < 1 (95 % 
C.L.) in Song et al. (2007a)) and are similar to the limits ob¬ 
tained in Marchini & Salvatelli (2013); Hu et al. (2013) (see 
also Dossett et al. (2014) where data from WiggleZ were used, 
Cataneo et al. (2015) where they considered galaxy clusters). 

5.3.4. Non-universal couplings: coupled Dark Energy 

Universal couplings discussed in the previous subsection gener¬ 
ally require screening mechanisms to protect baryonic interac¬ 
tions in high density environments, where local measurements 
are tightly constraining (see e.g. Khoury 2010 and Vikram et al. 
(2014) for astrophysical constraints). An alternative way to pro¬ 
tect baryons is to allow for non-universal couplings, in which 


■- P/anck+lensing 

P/anck+lensing+BSH 

- P/anck+lensing+WL 

P/anck+lensing+BAO/RSD 
- P/anck+lensing+BAO/RSD+WL 



log(B 0 ) 

Fig. 20. Likelihood plots of the f(R) theory parameter. Bo (see 
Sect. 5.3.3). CMB lensing breaks the degeneracy between Bo and 
the optical depth, t, resulting in lower upper bounds. 


different species can interact with different strengths: baryons 
are assumed to be minimally coupled to gravity while other 
species (e.g., dark matter or neutrinos) may feel a “fifth force,” 
with a range at cosmological scales. 

A fifth force between dark matter particles, mediated by the 
DE scalar field, is the key ingredient for the coupled DE sce¬ 
nario Amendola (2000). In the Einstein frame, the interaction is 
described by the Lagrangian 

-C = - V((f>) - m(tp + -CkintfA], (54) 

in which the mass of matter fields <// is not a constant (as 
in the standard cosmological model), but rather a function of 
the DE scalar field (f>. A coupling between matter and DE can 
be reformulated in terms of scalar-tensor theories or f(R ) mo¬ 
dels (Wetterich 1995b; Pettorino & Baccigalupi 2008; Wetterich 
2014) via a Weyl scaling from the Einstein frame (where matter 
is coupled and gravity is standard) to the Jordan frame (where the 
gravitational coupling to the Ricci scalar is modified and mat¬ 
ter is uncoupled). This is exactly true when the contribution of 
baryons is neglected. 

Dark matter (indicated with the subscript c) and DE densities 
are then not conserved separately, but coupled to each other: 

P0 = -yHp^X + w$) +f3p c cf)', (55) 

p' = -3Plp c - Pp c (/>'. 

Here each component is treated as a fluid with stress energy ten¬ 
sor T\ a)pl = ifi a + p a )u^u v + p a 6 v , where u^ = (-a, 0,0,0) is the 
fluid 4-velocity and w a = p a lp a is the equation of state. Primes 
denote derivatives with respect to conformal time and ji is as¬ 
sumed, for simplicity, to be a constant. This choice corresponds 
to a Lagrangian in which dark matter fields have an exponen¬ 
tial mass dependence ni((/>) = /;; () exp ^ (originally motivated 
by Weyl scaling scalar-tensor theories), where mo is a constant. 
The DE scalar field (expressed in units of the reduced Planck 
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f(R ) models 

Planck TT+lowP 

Planck TT+lowP 
+BSH 

Planck TT+lowP 
+WL 

Planck TT+lowP 
+BAO/RSD 

Planck TT+lowP 
+WL+BAO/RSD 

Bo . 

Bq (+lensing) . . . 

< 0.79 (95 % CL) 
<0.12 (95% CL) 

< 0.69 (95 % CL) 

< 0.07 (95 % CL) 

<0.10 (95% CL) 

< 0.04 (95 % CL) 

< 0.90 X 10~ 4 (95 % CL) 

< 0.97 x 10^ 4 (95 % CL) 

< 0.86 x 10^ 4 (95 % CL) 

< 0.79 x 10+95 % CL) 


Table 8 . 95 % CL intervals for the f(R) parameter, Bo (see Sect. 5.3.3). While the plots are produced for log I0 B ( >, the numbers in 
this table are produced via an analysis on Bo since the GR best fit value (By = 0) lies out of the bounds in a log 1() By analysis and its 
estimate would be prior dependent. 


mass M = (87rG ; y) l/2 j evolves according to the Klein-Gordon 
equation, which now includes an extra term that depends on the 
density of cold dark matter: 

4>" + I'Hcf)' + a 2 ^~ = a 2 /3p c . (56) 

d(f> 

Following Pettorino & Baccigalupi (2008), we choose an in¬ 
verse power-law potential defined as V — Vo<f> ", with a and 
Vo being constants. The amplitude Vo is fixed thanks to an 
iterative routine, as implemented by Amendola et al. (2012b); 
Pettorino et al. (2012). To a first approximation a only affects 
late-time cosmology. For numerical reasons, the iterative routine 
finds the initial value of the scalar field, in the range a > 0.03, 
which is close to the ACDM value a — 0 and extends the range 
of validity with respect to past attempts; the equation of state w 
is approximately related to a via the expression (Amendola et al. 
2012b): w = -2 /(a + 2) so that a value of a = 0.03 corresponds 
approximately to w(a = 0.03) = -0.99. The equation of state 
w = pip is not an independent parameter within coupled DE 
theories, being degenerate with the flatness of the potential. Dark 
matter particles then feel a fifth force with an effective gravita¬ 
tional constant G e ff that is stronger than the Newtonian one by a 
factor of /? 2 , i.e. 

G eff = G(l+2/? 2 ), (57) 

so that a value of/3 = 0 recovers the standard gravitational in¬ 
teraction. The coupling affects the dynamics of the gravitational 
potential (and therefore the late ISW effect), hence the shape 
and amplitude of perturbation growth, and shifts the position 
of the acoustic peaks to larger multipoles, due to an increase 
in the distance to the last-scattering surface; furthermore, it re¬ 
duces the ratio of baryons to dark matter at early times with 
respect to its present value, since coupled dark matter dilutes 
faster than in an uncoupled model. The strength of the coupling 
is known to be degenerate with a combination of G c , n s and ll(z) 
(Amendola & Quercellini 2003; Pettorino & Baccigalupi 2008; 
Bean et al. 2008; Amendola et al. 2012b). Several analyses have 
previously been carried out, with hints of coupling different from 
zero, e.g., by Pettorino (2013), who found f5 = 0.036 + 0.016 (us¬ 
ing Planck 2013 + WMAP polarization + BAO) different from 
zero at 2.2 <x (the significance increasing to 3.6 cr when data from 
FIST were included). 

The marginalized posterior distribution, using Planck 2015 
data, for the coupling parameter /? is shown in Fig. 21, while 
the corresponding mean values are shown in Table 9. Planck TT 
data alone gives constraints compatible with zero coupling and 
the slope of the potential is consistent with a cosmological con¬ 
stant value of cr = 0 at 1.3 cr. When other data sets are added, 
however, both the value of the coupling and the slope of the po¬ 
tential are pushed to non-zero values, i.e., further from ACDM. 
In particular, Planck+ BSH gives a value which is ~ 2.5 cr in 
tension with ACDM, while, separately, Planck+ WL+BAO/RSD 


gives a value of the coupling p compatible with the one from 
Planck +BSH and about 2.3 cr away from ACDM. 

When comparing with ACDM, however, the goodness of fit 
does not improve, despite the additional parameters. Only the 
Xbaoirsd i m P r °ves by ~ 1 in CDE with respect to ACDM, the 
difference not being significant enough to justify the additional 
parameters. The fact that the marginalized likelihood does not 
improve, despite the apparent 2 cr tension, may hint at some de¬ 
pendence on priors: for example, the first panel in Fig. 22 shows 
that there is some degeneracy between the coupling (5 and the 
potential slope a; while contours are almost compatible with 
ACDM in the 2 dimensional plot, the marginalization over a 
takes more contributions from higher values of /?, due to the 
degeneracy, and seems to give a slight more significant peak 
in the one dimensional posterior distribution shown in Fig. 21. 
Whether other priors also contribute to the peak remains to be 
understood. In any case, the goodness of fit does not point to¬ 
wards a preference for non-zero coupling. Degeneracy between 
the coupling and other cosmological parameters is shown in 
the other panels of the same figure, with results compatible 
with those discussed in Amendola et al. (2012b) and Pettorino 
(2013). Looking at the conservation equations (i.e., Eqs. 55 and 
56), larger positive values of fi correspond to a larger transfer 
of energy from dark matter to DE (effectively adding more DE 
in the recent past, with roughly Q,/, K P 1 for an inverse power- 
law potential) and therefore lead to a smaller Q m today; as a 
consequence, the distance to the last-scattering surface and the 
expansion rate are modified, with H' /H = -3/2(1 + w e ff), where 
w e ff is the effective equation of state given by the ratio of the to¬ 
tal pressure over total (weighted) energy density of the coupled 
fluid; a larger coupling prefers larger Hq and higher erg. 

The addition of polarization tightens the bounds on the cou¬ 
pling, increasing the tension with ACDM, reaching 2.8 cr and 
2.7 cr for Planck +BSH and Planck+ WL+BAO/RSD, respecti¬ 
vely. Also in this case the overall x~ does not improve between 
coupled DE and ACDM. 

6. Conclusions 

The quest for Dark Energy and Modified Gravity is far from 
over. A variety of different theoretical scenarios have been pro¬ 
posed in literature and need to be carefully compared with the 
data. This effort is still in its early stages, given the variety of the¬ 
ories and parameterizations that have been suggested, together 
with a lack of well tested numerical codes that allow us to make 
detailed predictions for the desired range of parameters. In this 
paper, we have provided a systematic analysis covering a gene¬ 
ral survey of a variety of theoretical models, including the use 
of different numerical codes and observational data sets. Even 
though most of the weight in the Planck data lies at high redshift, 
Planck can still provide tight constraints on DE and MG, espe¬ 
cially when used in combination with other probes. Our focus 
has been on the scales where linear theory is applicable, since 
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CDE models Planck TT+lowP Planck TT+lowP Planck TT+lowP Planck TT+lowP Planck TT+lowP 

+BSH +WL +BAO/RSD +WL+BAO/RSD 

p . <0.066(95%) 0.037!°-°]* 0.043!°“° 0.034!°°]]! 0.037!°“° 

0 43 + 015 n 7 q+°, 0 77 044 + 0.18 n 40+°- 15 n 45 +°- 17 

“. - 0.33 - 0.26 - 0.29 - 0.29 - 0.33 

Ho (km/s/Mpc) . . . . 65.4+3-2 67.47!°** 67.6 ± 2.8 66.7 ±1.1 66.9 ±1.1 

cr 8 . 0.812+°“! 0.829 ±0.018 0.819+°“! 0.817 ±0.017 0.810 ±0.017 


TT,TE,EE+lowP TT,TE,EE±lowP 
+BSH 


TT,TE,EE+lowP 

+WL 


TT,TE,EE+lowP TT,TE,EE+lowP 
+BAO/RSD ±WL±BAO/RSD 


P . <0.062(95%) 0.036+°°]° 

a . 0.42+°- 1 4 2 < 0.58 (95 %) 

Ho (km/s/Mpc) .... 65.4+i* 67.39+° *° 

o-g . 0.814+°“° 0.832 ±0.016 


0 036 +0019 
u.ujUq 026 

0 034+° ° 18 
u.uj+ 0,015 

0.038!°°]* 

"42-°;i; 

0 37+ 013 
- 0.28 

°- 42 -031 

66-5+2* 

66 7+' 1 

- 0.95 

66 - 8 -094 

0 81 7 + 0.031 
- 0.025 

0.823 ±0.015 

0.818 ±0.015 


Table 9. Marginalized mean values and 68 % C.L. intervals for coupled DE (see Sect. 5.3.4). Planck here refers to Planck TT+lowP. 
Results and goodness of fits are discussed in the text. CMB lensing does not improve the constraints significantly. 


Planck +BSH P!anck+\NL Planck+B AO/RSD Hi P/anc/t+WL+BAO/RSD 



0.2 0.4 0.6 0.8 0.25 0.30 0.35 0.40 0.76 0.80 0.84 0.88 63 66 69 72 


a <78 Hq 


Fig. 22. Marginalized posterior distribution for coupled DE and different combinations of the data sets (see Sect. 5.3.4). Here Planck 
refers to Planck TT+lowP. We show the degeneracy of the coupling f5 with a , Q m , erg and Hq. 



P 


Fig. 21. Marginalized posterior distribution for the coupling fi 
(see Sect. 5.3.4). The value corresponding to standard gravity is 
zero. Results and goodness of fit are discussed in the text. 


these are the most theoretically robust. Overall, the constraints 
that we find are consistent with the simplest scenario, ACDM, 
with constraints on DE models (including minimally-coupled 
scalar field models or evolving equation of state models) and 


MG models (including effective field theory, phenomenological 
parameterizations, f(R) and coupled DE models) that are signifi¬ 
cantly improved with respect to past analyses. We discuss here 
our main results, drawing our conclusions for each of them and 
summarizing the story-line we have followed in this paper to 
discuss DE and MG. 

Our journey started from distinguishing background and per¬ 
turbation parameterizations. In the first case, the background is 
modified (which in turn affects the perturbations), leading to the 
following results. 

1. The equation of state w(z) as a function of redshift has been 

tested for a variety of parameterizations. 

(a) In (wo,w a ), Planck TT+lowP+BSH is compatible with 
ACDM, as well as BAO/RSD. When adding WL to 
Planck TT+lowP, both WL and CMB prefer the (wo, w a ) 
model with respect to ACDM at about 2cr, although with 
a preference for high values of Hq (third panel of Fig. 3) 
that are excluded when including BSH. 

(b) We have reconstructed the equation of state in red- 
shift, testing a Taylor expansion up to the third orda¬ 
in the scale factor and by doing a Principal Component 
Analysis of w(z) in different redshifts bins. In addition, 
we have tested an alternative parametrization, that allows 
to have a varying w(z ) that depends on one parameter 
only. All tests on time varying w{z) are compatible with 
ACDM for all data sets tested. 
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2. ‘Background’ Dark Energy models are generally of 
quintessence type where a scalar field rolls down a poten¬ 
tial. We have shown via the (e s , e m ) parameterization, re¬ 
lated respectively to late and early time evolution, that the 
quintessence/phantom potential at low redshift must be rel¬ 
atively flat: d In V/dd) < 0.35 /Mp for quintessence; and 
d In V/cUb < 0.68/Mp for phantom models. A zero slope 
(ACDM) remains consistent with the data and compared to 
previous studies, the uncertainty has been reduced by about 
10 %. We have produced a new plot (Fig. 9) that helps to vi¬ 
sualize minimally coupled scalar field models, similarly to 
analogous plots often used in inflationary theories. 

3. Information on DE, complementary to (wo, vv a ), comes from 
asking whether there can be any DE at early times. First, 
we have obtained constraints on early DE parametrizations, 
assuming a constant DE relative density at all epochs until it 
matches ACDM in recent times: we have improved previous 
constraints by a factor ~ 3-4, leading to O c < 0.0071 at 95% 
C.L. from Planck TT+lowP+lensing+BSH and £2 e < 0.0036 
for Planck TT,TE,EE+lowP+BSH. In addition, we have also 
asked how much such tight constraints are weakened when 
the fraction of early DE is only present in a limited range in 
redshift and presented a plot of O e (z) as a function of z e , the 
redshift starting from which a fraction Q e is present. Also 
in this case constraints are very tight, with O e ^ 2 % (95 % 
C.L.) even for z e as late as ~ 50. 

The background is then forced to be very close to ACDM, unless 
the tight constraints on early DE can somehow be evaded in a 
realistic model by counter balancing effects. 

In the second part of the paper, we then moved on to un¬ 
derstanding what Planck can say when the evolution of the per¬ 
turbations is modified independently of the background, as is the 
case in most MG models. For that, we followed two complemen¬ 
tary approaches: one (top-down) that starts from a very general 
EFT action for DE (Sect. 5.2.1); the other (bottom-up) that starts 
from parameterizing directly observables (Sect. 5.2.2). In both 
cases we have assumed that the background is exactly ACDM, 
in order to disentangle the effect of perturbations. We summarize 
here our results. 

1. Starting from EFT theories for DE, which include (almost) 
all universally coupled models in MG via nine generic fun¬ 
ctions of time, we have discussed how to restrict them to 
Horndeski theories, described in terms of five free functions 
of time. Using the publicly available code EFTcamb, we have 
then varied three of these functions (in the limits allowed 
by the code) which correspond to a non-minimally coupled 
K-essence model (i.e. a M , a K are varying functions of 
the scale factor). We have found limits on the present value 
Q'mo < 0.052 at 95% C.L. (in the linear EFT approximation), 
in agreement with ACDM. Constraints depend on the stabi¬ 
lity routines included in the code, which will need to be fur¬ 
ther tested in the future, together with allowing for a larger 
set of choices for the Horndeski functions, not available in 
the present version of the numerical code. 

2. When starting from observables, two functions of time and 
scale are required to describe perturbations completely, in 
any model. Among the choices available (summarized in 
Sect. 3.2.2), we choose fi(a,k ) and rj(a,k) (other observ¬ 
ables can be derived from them). In principle, constraints 
on these functions are dependent on the chosen parame¬ 
terization, which needs to be fixed. We have tested two 
different time dependent parameterizations (DE-related and 
time-related) and both lead to similar results, although the 


first is slightly more in tension than the other with ACDM 
(Fig. 14). In this framework, ACDM lies at the 2 cr limit when 
Planck TT+lowP+BSH is considered, the tension increa¬ 
sing to about 3cr when adding WL and BAO/RSD to Planck 
TT+lowP. As discussed in the text, the mild tension with 
Planck TT+lowP is related to lower power in the TT spec¬ 
trum and a larger lensing potential in the MG model, with 
respect to ACDM. The inclusion of CMB lensing shifts all 
contours back to ACDM. We have reconstructed the two ob¬ 
servables in redshift for both parameterizations, along the 
maximum degeneracy line (Fig. 16). When scale dependence 
is also included, constraints become much weaker and the 
goodness of fit does not improve, indicating that the data do 
not seem to need the addition of additional scale dependent 
parameters. 

The last part of the paper discusses a selection of particular MG 
models of interest in literature. 

1. We first commented on the simple case of a minimally cou¬ 
pled scalar field in which not only the equation of state is 
allowed to vary but the sound speed of the DE fluid is not 
forced to be 1, as it would be in the case of quintessence. 
Such a scenario corresponds to k-essence type models. As 
expected, given that the equation of state is very close to the 
ACDM value, the total impact of DE perturbations on the 
clustering is small. 

2. We adopt an alternative way to parameterize observables 
(the equation of state approach) in terms of gauge invari¬ 
ant quantities T and cr. We have used this approach to in¬ 
vestigate Lorentz-breaking massive gravity and generalised 
scalar fields models, updating previous bounds. 

3. As a concrete example of universally coupled theories, we 
have considered f(R) models, written in terms of B(z), con¬ 
ventionally related to the first and second derivatives of f(R) 
with respect to R. Results are compatible with ACDM. Such 
theories assume that some screening mechanism is in place, 
in order to satisfy current bounds on baryonic physics at so¬ 
lar system scales. 

4. Alternatively to screening mechanisms, one can assume that 
the coupling is not universal, such that baryons are still feel¬ 
ing standard gravitational attraction. As an example of this 
scenario, we have considered the case in which the dark 
matter evolution is coupled to the DE scalar field, feel¬ 
ing an effective fifth force stronger than gravity by a fac¬ 
tor /j 1 2 . Constraints on coupled dark energy show a tension 
with ACDM at the level of about 2.5 cr, slightly increasing 
when including polarization. The apparent tension, however, 
seems to hint at a dependence on priors, partly related to the 
degeneracy between the coupling and the slope of the back¬ 
ground potential (and possibly others not identified here). 
Future studies will need to identify the source of tension and 
possibly disentangle background from perturbation effects. 

There are several ways in which the analysis can be ex¬ 
tended. We have made an effort to (at least start) to put some 
order in the variety of theoretical frameworks discussed in lit¬ 
erature. There are of course scenarios not included in this pic¬ 
ture that deserve future attention, such as additional cosmolo¬ 
gies within the EFT (and Horndeski) framework, Galileons (see 
for example Barreira et al. (2014)), other massive gravity models 
(see de Rham (2014) for a recent review), general violations of 
Lorentz invariance as a way to modify GR (Audren et al. 2014), 
non-local gravity (which, for some choices of the action, ap¬ 
pears to fit Planck 2013 data sets (Dirian et al. 2014) as well as 
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ACDM, although there is no connection to a fundamental the¬ 
ory available at present); models of bigravity (Hassan & Rosen 
2012) appear to be affected by instabilities in the gravitational 
wave sector (Cusin et al. 2014) (see also Akrami et al. (2015)) 
and are not considered in this paper. In addition to extending the 
range of theories, which requires new numerical codes, future 
tests should verify whether all the assumptions (such as stabil¬ 
ity constraints, as pointed out in the text) in the currently avail¬ 
able codes are justified. Further promising input may come from 
data sets such as WL and BAO/RSD, that allow to tighten con¬ 
siderably constraints on MG models in which perturbations are 
modified. With the data available at the time of this paper, there 
seems to be no significant trend (at more than 3 standard de¬ 
viations) that compensates any possible tension in erg or Hq by 
favouring Modified Gravity; nevertheless, this issue will need to 
be further investigated in the future. We also anticipate that these 
constraints will strengthen with future releases of the Planck 
data, including improved likelihoods for polarization and new 
likelihoods, not available at the time of this paper, such as ISW, 
ISW-lensing and fi-mode polarization, all of which can be used 
to further test MG scenarios. 
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